A NEW  SEMIEMPIRICAL  METHOD  FOR  ELECTRONIC  STRUCTURE 


By 

YA-WEN  HSIAO 


A DISSERTATION  PRESENTED  TO  THE  GRADUATE  SCHOOL 
OF  THE  UNIVERSITY  OF  FLORIDA  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE  DEGREE  OF 
DOCTOR  OF  PHILOSOPHY 


UNIVERSITY  OF  FLORIDA 


1999 


To  my  dear 
father, 

little  sister,  Jing-Wen, 
husband,  Magnus, 

and  to  the  memory  of  my  dear  mother 


ACKNOWLEDGMENTS 


I would  like  to  thank  my  advisor,  Professor  Michael  C.  Zerner,  for  bringing  this 
challenging  project  to  my  attention,  providing  great  help  during  my  years  of  study, 
and  many  valuable  suggestions  to  my  thesis. 

Many  thanks  to  the  members  of  my  committee-Professors  N.  Yngve  Ohrn,  John 
R.  Sabin,  Philip  J.  Brucat,  and  Nigel  G.  J.  Richards-for  their  comments  on  my  thesis, 
which  have  been  most  helpful. 

Special  thanks  to  Professor  Hendrik  J.  Monkhorst  for  being  willing  to  stand  in 
as  examiner  during  both  my  oral  and  final  defenses. 

I greatly  appreciate  the  daily  advice  and  friendship  of  Drs.  Marshall  Cory, 
Junqiang  Sun,  and  Ajith  Perera,  since  I joined  QTP. 

Immense  appreciation  is  dedicated  to  my  husband,  Magnus,  for  all  his  support, 
discussion,  and  love.  Also  I thank  my  father  and  sister  for  their  constant  encourage- 
ment, understanding,  and  comfort.  Without  them  I simply  cannot  finish  anything. 


TABLE  OF  CONTENTS 


page 

ACKNOWLEDGMENTS iii 

ABSTRACT  vii 

CHAPTERS 

1 INTRODUCTION 1 

2 THEORY  REVIEW  5 

2.1  The  Molecular  Problem 5 

2.1.1  Time-independent  Schrbdinger  Equation 5 

2.1.2  Born-Oppenheimer  Approximation 6 

2.1.3  Variational  Principle 7 

2.1.4  The  Hartree-Fock  Method  for  Molecules 8 

2.1.5  Koopmans’  Theorem 13 

2.1.6  Electron  Correlation 15 

2.2  Semiempirical  Methods 16 

2.2.1  Hiickel  and  Extended  Hiickel  Methods 17 

2.2.2  Semiempirical  SCF  MO  Methods  of  the  ZDO  Type 20 

2.2.3  Parameterization  for  the  INDO  Method  30 

2.2.4  General  Remarks 41 

3 XINDO 44 

3.1  Features 44 

3.2  Details  of  Building  XINDO 46 

3.2.1  Comparison  with  Previous  Models 46 

3.2.2  Integral  by  Integral 49 

3.2.3  Details  of  the  Parameterization 63 

3.2.4  Similar  Work  in  the  Field  71 

4 PARAMETERS  AND  OTHER  RESULTS 76 


IV 


page 


4.1  Parameters  7g 

4.2  Test  Cases 81 

4.2.1  Geometry  Optimization 81 

4.2.2  C3H3+  ; ; ; 82 

4.2.3  O3 83 

5 CONCLUDING  REMARK  85 

REFERENCES 88 

BIOGRAPHICAL  SKETCH 92 


V 


Abstract  of  Dissertation 
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A NEW  SEMIEMPIRICAL  METHOD  FOR  ELECTRONIC  STRUCTURE 

By 

Ya-Wen  Hsiao 
1999 


Chairman:  Michael  C.  Zerner 
Major  Department:  Chemistry 

A modified  semiempirical  molecular  orbital  method,  XINDO,  is  developed  at 
the  INDO  level.  The  commonly  applied  Goeppert-Mayer-Sklar  approximation  is  no 
longer  included,  and  a short  range  repulsive  term,  the  effective  core  potential,  is 
introduced  to  compensate  the  penetration  effect  which  causes  errors  in  calculating 
molecular  geometries.  The  XINDO  Hamiltonian  is  constructed  over  the  Slater  type 
functions.  During  each  SCF  cycle,  we  solve  the  secular  equation  with  non-identity 
overlap  matrix.  In  this  way,  we  incorporate  the  overlap  into  the  XINDO  method.  This 
improves  on  the  geometry  calculations,  and  the  energy  description  of  antibonding 
molecular  orbitals.  Multi-center  terms  are  also  added  by  both  modifying  the  off- 
diagonal  part  of  the  STO  Hamiltonian  representation,  and  indirectly  through  the 
SCF  procedure. 


VI 


The  XINDO  model  has  been  implemented  for  H,  C,  N,  and  O which  are  the 
main  constituents  of  most  of  the  organic  molecules.  The  model  can  be  applied  to 
calculate  electronic  excitations  and  molecular  geometries  by  using  the  semiempirical 
two-electron  integrals  which  are  conventionally  used  only  for  calculating  spectra.  The 
results  of  the  calculated  geometries  are  quite  satisfactory,  and  the  calculated  charge 
distributions  and  inter-atomic  angles  are  improved  through  the  inclusion  of  multi- 
center integrals. 


CHAPTER  1 
INTRODUCTION 


After  several  decades,  the  development  of  electronic  structure  computational 
techniques  has  become  quite  mature.  Ab  initio  methods  provide  accurate  results  for 
smaller  molecular  systems,  especially  for  ground  state  properties.  However,  when  it 
comes  to  systems  of  biological  or  pharmacological  significance,  ab  initio  treatment 
becomes  prohibitive  due  to  the  computational  demands  that  arise  from  the  increase 
in  the  number  of  electrons.  Moreover,  a reasonable  description  of  the  excited  states 
of  molecules  is  still  difficult  to  achieve.  However,  it  is  possible  today  to  use  ab  initio 
methods  to  do  expensive  calculations  to  obtain  approximate  excited  state  descriptions 
for  small-sized  molecules.  Nevertheless,  ab  initio  methods  are  still  not  applicable  in 
practical  studies  involving  large  systems.  For  example,  the  study  of  the  primary 
electron  transfer  mechanism  of  bacterial  photosynthesis  is  an  active  subject.  The 
electron  moves  from  the  special  pair  to  a bacteriopheophytin-6  within  3 ps  and  hence 
makes  experimental  investigations  difficult.  Since  the  entire  protein  complex  consists 
of  more  than  10,000  heavy  atoms,  it  is  impractical  to  use  ab  initio  quantum  mechanics 
for  theoretical  predictions. 
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Semiempirical  electronic  structure  methods,  on  the  other  hand,  offer  an  alter- 
native way  of  treating  the  electronic  structures  of  large  systems  using  quantum  me- 
chanics. Most  of  the  large  proteins  or  complex  compounds  consist  of  smaller-sized 
functional  groups  in  which  individual  molecular  properties  are  conserved.  Based  on 
these  smaller  organic  molecules  many  semiempirical  methods  are  parameterized  to  fit 
the  experimental  data  of  the  these  properties.  This  offers  an  advantage  to  semiempir- 
ical methods  in  that  they  are  able  to  describe  all  classes  of  molecules  averagely  well 
for  general  purposes.  The  disadvantage  of  using  the  semiempirical  methods  is  that 
calculations  of  electronic  spectra  are  not  available  in  most  of  the  methods,  and  the 
methods  which  are  able  to  calculate  spectra  cannot  in  general  be  used  to  calculate 
molecular  geometry.  Therefore,  the  few  methods  for  calculating  electronic  excitations 
do  not  give  meaningful  ground  state  potential  energy  surfaces,  not  to  mention  excited 
states  potential  energy  surfaces. 

The  motivation  of  this  thesis  is  to  find  a semiempirical  method  which  can  yield 
reasonable  molecular  geometries  as  well  as  electronic  spectra. 

There  may  arise  two  questions  to  this  motivation.  First  of  all,  why  do  we  choose 
semiempirical  methods?  Because  we  want  to  have  a method  which  can  help  in  our 
studies  on  big  molecular  systems,  and  as  it  has  been  stated  previously,  ab  initio 
methods  are  still  restricted  to  smaller  molecules.  The  second  question  is:  why  do 
we  bother  to  have  a method  which  is  multi-purposed?  The  most  straightforward 
answer  is  that  such  a method  is  more  convenient  for  users  to  use  for  calculations  of 
all  properties,  instead  of  changing  from  method  to  method  only  because  the  range 
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of  applicability  of  each  methods  is  limited.  Another  advantage  is  that  the  calculated 
properties  are  consistent.  The  ground  state  potential  energy  surface  obtained  from 
such  method  will  be  correct  within  the  accuracy  of  the  method,  and  the  electronic 
excitations  will  also  be  good.  Then  the  next  possible  step  is  to  create  the  excited 
state  potential  energy  surface  based  on  the  ground  state  potential  surface  and  the 
electronic  excitations.  With  such  a method  one  may  imagine  that  it  becomes  possible 
to  study  the  excited  states  of  a system  as,  e.g.,  the  photosynthesis  reaction  center 
mentioned  previously. 

Yet,  there  are  other  problems  existing  in  semiempirical  methods,  such  as  the 
difficulty  in  obtaining  correct  conformational  structures,  rotational  barriers,  binding 
energies,  etc..  This  makes  the  semiempirical  methods  open  to  further  improvement. 
Hence  we  will  attempt  to  design  the  new  method  to  improve  these  shortcomings  as 
they  become  known. 

Many  assumptions  have  been  made  for  semiempirical  methods  in  order  for  them 
to  yield  fast  and  correct  answers  with  low  cost.  In  the  new  method,  however,  we  will 
bring  back  more  of  the  physics  that  has  been  omitted  in  earlier  semiempirical  meth- 
ods. More  integrals  may  need  to  be  added,  but  this  must  be  done  with  great  care  so 
that  we  will  not  destroy  the  advantage  of  the  semiempirical  methods.  We  first  review 
the  most  commonly  used  semiempirical  methods:  the  ZDO  (zero  differential  overlap) 
type  methods  and  the  EHT  (extended  Hiickel  theory).  Especially  the  assumptions 
behind  the  methods  need  to  be  understood  and  the  invalidity  re-examined.  Second, 
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we  will  learn  from  the  existing  models  and  combine  the  useful  parts  of  them  into  the 
new  model. 

In  summary,  we  think  that  we  have  adopted  the  strong  points  of  both  the 
EHT  method  (the  inclusion  of  the  overlap)  and  the  ZDO  methods  (systematically 
reducing  the  number  of  two-electron  integrals)  in  the  design  of  XINDO.  The  molecular 
geometries  obtained  by  using  the  XINDO  method  are  quite  satisfactory.  The  choice 
of  the  two-electron  integrals  promises  the  applicability  of  the  electronic  excitation 
calculations.  The  new  method,  XINDO,  is  indeed  full  of  potential. 


CHAPTER  2 
THEORY  REVIEW 

2.1  The  Molecular  Problem 
2.1.1  Time-independent  Schrodinger  Equation 

Solving  the  time-independent  Schrddinger  equation  [1] 


= Ej^i 


(2.1) 


offers  a way  to  obtain  properties  of  chemical  substances.  In  the  Schrodinger  equation, 
H is  the  Hamiltonian  operator,  Ei  is  the  total  energy  of  the  state  I of  the  system, 
and  is  the  wavefunction  describing  the  stationary  states  of  the  system,  e.g.,  a 
molecule.  The  usual  time-independent  Hamiltonian  operator  for  a molecular  system 
consisting  of  N electrons  and  M atoms  is 


H = 
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A-KEaRAB 


(2.2) 


where  the  first  two  sums  are  the  kinetic  energy  operators  of  electrons  and  nuclei, 
respectively,  the  third  sum  is  the  attraction  between  electrons  and  nuclei,  and  the  last 
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two  sums  are  the  electron-electron  and  nuclear-nuclear  repulsions,  respectively.  Other 
interactions  such  as  spin-orbit  coupling  can  be  added,  depending  on  the  purpose  of 
the  calculation. 

2.1.2  Born-Oppenheimer  Approximation 

Except  for  the  nuclear-electron  attraction  term  in  Eqn.  (2.2),  other  electronic 
and  nuclear  terms  are  independent  of  each  other.  It  seems  possible  to  simplify  the 
process  of  solving  the  Schrodinger  equation  by  separating  the  electronic  and  nuclear 
motions,  i.e.,  the  Born-Oppenheimer  [2]  approximation.  This  approximation  is  based 
on  the  fact  that  the  nuclear  mass  is  at  least  1836  times  the  mass  of  an  electron, 
hence  the  electrons  move  much  faster  than  the  nuclei.  One  may  therefore  assume 
that  electrons  move  in  a field  generated  by  a frozen  nuclear  configuration.  The 
approximate  molecular  wavefunction  is  written  as  a product  of  an  electronic  part 
and  a nuclear  part 

^(x,  .Y)  = X)'il)nud{X),  (2.3) 

where  x indicates  the  electronic  coordinates  and  is  the  nuclear  coordinates.  The 
notation  tl;ei{x-,X)  indicates  that  the  nuclear  coordinates  enters  the  electronic  wave- 
function  only  parametrically.  For  a given  nuclear  configuration.  A",  the  electronic 
wavefunction  may  then  be  obtained  from  solving  the  electronic  part  of  the  Schrbdinger 
equation. 


Helix;  X)^P,i{x;X)  = Eet{X)M^;X). 


(2.4) 
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Eei{X)  together  with  the  nuclear  repulsion  energy,  which  is  a constant  for  each  set  of 
and  can  therefore  be  added  after  Eqn.  (2.4)  is  solved,  serve  as  a potential  energy 
Kiud(A^)  in  which  the  nuclei  move  [2,3].  The  global  minimum  of  Vnud{X)  occurs  at 
the  most  stable  nuclear  configuration,  i.e.,  at  the  equilibrium. 

2.1.3  Variational  Principle 

Consider  a system  described  by  a Hamiltonian  H and  with  a normalized  ground 
state  wavefunction  Then  the  energy  expectation  value  of  the  Hamiltonian  with 
respect  to  any  normalized  trial  wavefunction  $,  is  always  greater  or  equal  to  the  true 
ground  state  energy  of  the  system: 

= ($|//|$)  > Eo.  (2.5) 

The  energy  £'($)  is  a functional  of  the  $’s,  and  fully  minimizing  £"($)  with  respect 
to  variations  in  the  trial  wavefunctions  gives  the  true  ground  state  and  the  energy 
£’(^o)  = Eo,  which  is  indicated  in  Eqn.  (2.5)  where  the  equality  holds  when  $ is 
identical  to  ^o- 

Since  the  Schrodinger  equation  cannot  be  solved  exactly,  except  for  the  simplest 
cases,  we  are  more  interested  in  finding  approximate  solutions.  One  way  is  to  apply 
the  variational  principle,  since  one  may  always  find  an  upper  bound  to  Eq  according 
to  Eqn.  (2.5). 
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2.1.4  The  Hartree-Fock  Method  for  Molecules 

Due  to  the  cumbersome  electron-electron  interaction  term,  the  only  systems  of 
electrons  and  nuclei  that  can  be  exactly  solved  are  the  hydrogen  atom  and  systems 
such  as  within  the  Born-Oppenheimer  approximation  because  there  is  only  one 
electron.  For  many-electron  systems,  since  the  inter-electronic  repulsion  term  in  the 
Schrodinger  equation  is  not  separable,  further  approximation  is  needed.  One  can 
obtain  the  zeroth  order  wavefunction  of  any  many-electron  system  by  neglecting  the 
inter-electronic  repulsions,  and  then  the  solution  to  the  resulting  Schrbdinger  equation 
can  be  expressed  as  a product  of  many  one-electron  functions: 

^(xi, X2,  • • • , X;\f)  = ipi{ri  : Si)ip2{T^2  ■ S2)  ■ • • (2.6) 

where  Xi  contains  the  spatial  (r,)  and  spin  (sj)  descriptions  of  the  ith  electron,  i.e., 
Xi  = {ti  : Si).  The  procedure  of  calculating  these  one-electron  functions  was  intro- 
duced by  Hartree  [4].  In  the  Hartree  theory  the  product  wavefunction  Eqn.  (2.6) 
is  used  as  a variational  many-electron  wavefunction  in  Eqn.  (2.5).  This  results  in  a 
Schrodinger  like  equation  for  each  ipi  where  the  cumbersome  electron-electron  inter- 
action is  replaced  by  an  average  interaction  (same  as  Eqn.  (2.14)). 

From  now  on  we  will  proceed  our  discussion  using  atomic  units  in  which  h, 
e^j4TT£o  and  rUe  are  all  equal  to  one. 

However,  since  the  true  many-electron  Schrodinger  equation  is  not  separable,  the 
true  wavefunction  can  not  be  expressed  exactly  as  a product  of  many  one-electron 
functions  and,  moreover,  Eqn.  (2.6)  does  not  satisfy  the  Pauli  exclusion  principle 
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[5],  which  states  that  a many-fermion  wavefunction  has  to  change  sign  under  the 
interchange  of  two  particles, 


(^1 ) ‘ ‘ ) ' ' ' ) ‘ ‘ j ) 


(^1 ) ' ‘ ‘ ‘ ‘ 5 ' ■ ■ ? • 


(2.7) 


The  state  of  an  electron  is  not  determined  from  the  spatial  orbital  only  but  one 
also  has  to  consider  its  spin  state.  Keeping  the  single  particle  approximation,  the 
approximate  many-electron  wavefunction  must  be  written  as  an  antisymmetric  linear 
combination  of  product  wavefunctions  using  spin  orbitals.  This  was  pointed  out  by 
Fock  [6]  and  Slater  [7].  An  SCF  calculation  that  uses  antisymmetrized  spin  orbitals 
is  called  a Hartree-Fock  (HF)  calculation.  A convenient  way  of  representing  the 
antisymmetric  linear  combination  of  the  Hartree  SCF  wave  function,  is  the  Slater 
determinant  [8] 


^ = 


\/M 


t/'l(l)  V'2(l) 

v>i(2)  V’2(2) 

MN)  MN) 


V’iv(2) 

^n{N) 


(2.8) 


where  V are  the  permutation  operators  of  the  symmetric  group  of  order  N. 

With  the  assumption  that  the  Slater  determinant  Eqn.  (2.8)  is  an  approximation 
to  the  ground  state  many-electron  wavefunction,  one  may  use  the  variational  principle 
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Eqn.  (2.5)  to  minimize  the  energy  given  by 


E,  = E('l'*W  + 5Efe1ly> 

= i (2.9) 

* i,j 

where  the  one  electron  operator  h is  given  by 


= (2.10) 

2 ^ riA 

and  the  two-electron  integrals  are  given  by 

(ij\kl)  = j dxidx2ip*{l)'ip*{2)r^2^^pk{l)i;i{2).  (2.11) 

Such  a procedure  leads  to  a Schrodinger-like  equation  for  the  one-electron  wave 
functions,  the  Hartree-Fock  equation; 

M = Eii/Ji,  (2.12) 


where  the  Fock-operator  is  given  by 


/(I) 


3 

OCC  - 


(2.13) 
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The  operator  Vu  interchange  electron  1 and  2,  and  Jj  is  the  Coulomb  operator,  which 
represents  the  average  Coulomb  interaction  between  the  electrons.  Kj  is  the  exchange 
operator,  which  comes  from  the  requirement  of  an  anti-symmetric  wavefunction;  that 
is. 

Coulomb  operator: 


Jj^i{l)  = J da;2V’J(2)ri2Vt(l)^j(2),  (2.14) 

Exchange  operator: 

= j dx2t/'*(2)V^i(2)rf2Vj(l)  (2.15) 

From  Eqn.  (2.13)  it  is  clear  that  the  Fock  operator  depends  on  the  wavefunction 
through  the  Coulomb  and  exchange  operators.  Therefore  the  Hartree-Fock  equation 
is  not  linear  and  is  generally  solved  iteratively. 

A particular  way  of  solving  the  HF  equation  is  to  write  each  molecular  orbital  as 
a linear  combination  of  atomic  orbitals  (LCAO-MO)  [9, 10].  Let  {0i,  • • • , ■ • • , 

denote  a set  of  K real  atomic  orbital  functions  to  build  up  our  trial  zth  molecular 
orbital  function  %l)i 

K 

fi=l 


(2.16) 
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where  c^^’s  are  the  expansion  coefficients  to  be  determined.  Solving  the  Hartree-Fock 
equation  Eqn.  (2.12)  with  this  trial  wavefunction,  one  obtains 

/(i)v^,(i) = /(I)  f:  = Ei  x:  (2.17) 

/J=l  ^L=l 

Multiplying  Eqn.  (2.17)  from  the  left  by  (/>*  and  integrate,  one  arrives  at  the  secular 
equation,  often  called  Roothaan  equation, 

Fc  = AcE,  (2.18) 

where  E is  a diagonal  matrix  with  the  energy  eigenvalues  on  the  diagonal  and  A = 
1 + S.  The  matrices  S and  F are  defined  through 

= (0^(1)|0,.(1)),  for  /i  ^ z/,  (2.19) 

= (<^^(1)|/(1)|0^(1)).  (2.20) 

From  Eqn.  (2.18)  we  get  K eigenfunctions,  • • • , V’o  • • • > '0/c,  with  the  correspond- 
ing energy  eigenvalue  e,; 

occ 

Si  = {i\h\i)  + Y,{ij\\i3)-  (2.21) 


In  the  case  of  closed-shell  A^-electron  system,  N/2  spatial  eigenfunctions  are  occupied. 
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The  approximate  electron  energy  E^i, 


Eel  = 


(2.22) 


is  the  upper  bound  of  the  true  ground  state  energy,  according  to  the  variational 
principle.  The  energy  is  not  equal  to  the  sum  of  eigenvalues,  since  in  that  case 
the  two-electron  interaction  would  be  counted  twice.  As  was  mentioned  above,  the 
Fock  equation  depends  on  the  trial  functions,  and  is  generally  solved  by  an  iterative 
process.  This  is  called  the  self-consistent-field  method.  The  criteria  for  deciding 
when  to  stop  the  iteration  is  most  commonly  given  by  the  difference  of  the  density 
matrices  of  the  current  cycle  and  the  previous  one.  The  definition  of  the  density 
matrix  elements  is 


where  Ui  is  the  occupation  number  of  molecular  orbital  i.  If  the  density  matrix  is 
equal  to  the  result  from  the  previous  iteration  within  a given  tolerance,  then  the 
electrostatic  field  is  said  to  be  self-consistent. 

2.1.5  Koopmans’  Theorem 

In  the  previous  section  it  was  shown  that  adding  up  the  HF  orbital  energies  6i 
does  not  give  the  total  energy,  since  in  that  case  the  Coulomb  and  exchange  energies 
would  be  counted  twice.  What  physical  meaning  can  be  given  to  the  HF  orbital 


(2.23) 


energies? 
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Consider  an  iV-electron  system  with  a ground  state  Slater  determinant 
built  from  the  A^-HF  orbitals  ijji  with  the  lowest  energies.  Take  this  A^-electron 
system  and  compare  to  the  systems  where  an  electron  is  either  added  or  removed, 
but  keeping  the  set  of  orbitals  ipi,  - ■ • obtained  from  the  HF  calculation  on  the 
A^'-electron  system.  The  energy  expectation  values  are 


^Eo  = 

(2.24) 

(2.25) 

N+l^a  _ ^Af+l^a| 

(2.26) 

where  the  subscript  i denotes  that  an  electron  has  been  removed  from  occupied  orbital 
tpi  and  the  superscript  a means  that  an  electron  is  added  to  a virtual  orbital  0a-  The 
ionization  potential  in  this  frozen  orbital  frame: 

Eo  = (2.27) 

and  similarly, 

^Eq  -^+1  E^^  = -6a,  (2.28) 

is  the  electron  affinity. 

However,  Koopmans’  theorem  tends  to  give  too  positive  ionization  potentials 
and  too  negative  electron  affinities.  Yet,  Koopmans’  approximation  is  still  reasonable 
for  the  ionization  potential  for  a neutral  species  or  a positive  ion  for  delocalized 
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molecular  orbitals,  although  it  is  not  ideal  for  the  electron  affinity  of  a neutral  or  a 
negative  species.  There  are  two  error  sources  when  applying  Koopmans’  theorem. 
The  first  is  that  the  orbitals  of  the  ionic  system  are  assumed  to  be  the  same  as 
the  orbitals  of  the  neutral  system;  i.e.,  the  orbitals  are  frozen.  The  energy  for  the 
ion  calculated  with  these  orbitals  is  higher  than  it  would  be  if  one  would  allow  for 
the  orbitals  to  relax.  The  other  source  of  error  is  the  lack  of  correlation  in  the  HF 
procedure.  One  would  expect  the  correlation  correction  to  be  larger  for  the  state  with 
the  larger  number  of  electrons.  Since  the  correlation  energy  is  a negative  quantity 
one  may  assume  that  the  neglect  of  electron  correlation  somewhat  cancel  the  error 
from  the  frozen  orbitals  in  the  case  of  ionization  potentials  but  adds  to  the  error  in 
the  estimate  of  the  electron  affinities. 

2.1.6  Electron  Correlation 

A Hartree-Fock  SCF  wave  function  describes  the  interactions  among  electrons 
in  an  average  way  as  shown  in  Eqn.  (2.14)  . The  error  in  the  energies  from  HF  calcu- 
lations is  of  the  magnitude  of  1%  which  can  be  very  dramatic  for  certain  properties, 
such  as  bonding  energies. 

The  difference  between  the  exact  energy  for  a specified  non-relativistic  Hamil- 
tonian and  the  HF  energy  for  the  same  system  is  defined  as  the  correlation  energy 
[11]. 

What  is  the  nature  of  the  correlation  energy?  An  electron  in,  e.g.,  a molecule 
feels  the  actual  Coulomb  repulsion  from  the  other  electrons  and  therefore  tend  to 
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avoid  each  other  more  than  what  is  the  case  when  using  the  average  Coulomb  poten- 
tial Eqn.  (2.14).  In  other  words  the  electron  motion  is  correlated  and  this  gives  rise  to 
a lowering  of  the  energy.  This  is  normally  referred  as  the  “dynamic  correlation”  and 
it  usually  represents  the  largest  part  of  the  total  correlation  energy.  The  other  source 
of  the  electron  correlation  is  the  so-called  “essential  correlation.”  It  arises  when  one 
attempts  to  describe  the  electronic  structures  with  a single  Slater  determinant. 

Solving  the  full  Schrodinger  equation  with  the  actual  electron  repulsion  is  im- 
possible for  a system  with  many  electrons.  However,  there  are  other  ways  to  at  least 
partially  calculate  the  correlation  energy  and  thus  correct  the  HF  results.  Two  such 
methods  are  configuration  interaction  (Cl),  and  many-body  perturbation  theory. 

2.2  Semiempirical  Methods 

Despite  the  lack  of  electron  correlation,  the  HF  SCF  method  is  quite  accurate  in 
describing  molecular  properties.  However,  the  evaluation  of  integrals  has  made  the 
method  limited  to  smaller  molecules  and  highly  demanding  of  computer  resources. 
The  cost  of  calculation  increases  as  the  number  of  electrons  increases.  Given  K basis 
functions  then  the  evaluation  of  two-electron  integrals  is  formally  an  A''*  process,  and, 
in  order  to  get  better  description  of  the  system,  large  basis  sets  are  needed.  In  order 
to  shed  some  light  of  large  systems  like  proteins,  organo-metalic  complex  compounds, 
alternative  ways  must  be  sought. 
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Consider  the  molecular  Hamiltonian  expressed  in  second  quantization  formalism 
in  an  orthonormal  molecular  basis  {ipp}  [12]: 

^ + 7 (2.29) 

p,g  ^ pqrs 

where  a+  is  a creation  operator  which  creates  an  electron  in  spin  orbital  ^p,  and  a, 
is  an  annihilation  operator  which  removes  an  electron  from  spin  orbital  ipq.  These 
creation  and  annihilation  operators  fulfill  the  anticommutation  relations 

[®p  > ]+  0,  [dp,  dg]+  — 0,  [dp,  d^  = 5pq.  (2.30) 

By  observing  Eqn.(2.29)  and  one  can  see  that  the  integrals  hpq  and  {pq\\rs)  appear  as 
parameters.  If  the  choice  of  parameters  is  based  on  Eqs.(2.13),  (2.14),  and  (2.15)  then 
the  procedure  is  called  ab  initio.  On  the  other  hand,  one  can  find  suitable  parameters 
to  replace  the  integrals  in  the  ab  initio  method,  and  yet  give  results  of  the  molecular 
properties  which  fit  best  to  either  ab  initio  or  experimental  results.  Such  approaches 
are  called  semiempirical.  Since  many  integrals  can  be  replaced  with  parameters  in 
semiempirical  approaches,  it  is  quite  possible  to  do  calculations  on  large  molecules 
with  such  methods. 

2.2.1  Hfickel  and  Extended  Hfickel  Methods 

The  first  semiempirical  method  was  developed  by  Hfickel  [13]  in  1931  and  is 
referred  to  as  the  Hfickel  molecular-orbital  (HMO)  method.  HMO  consists  of  a one- 
electron  Hamiltonian,  that  takes  only  valence  7r-electrons  into  account.  It  differs  from 
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an  SCF  scheme,  since  the  electron  operators  do  not  depend  on  the  charge  density. 

= (2.31) 

i=l 

where  n is  the  number  of  tt  electrons.  The  form  of  is  not  specified,  but  it 

satisfies  the  relation 

= Eitpi,  (2.32) 

where  -ipi  is  the  zth  tt  molecular  orbital,  and  Si  is  the  molecular  orbital  energy.  HMO 
also  assumes  LCAO-MO  for  -tpi: 

M 

— 'y  (2.33) 

fl=l 

where  M is  the  number  of  atoms  and  is  a tt  atomic  orbital  on  the  /ith  atom.  Then 
the  variational  theorem  is  applied  and  the  molecular  orbitals  are  obtained  by  solving 
the  secular  equation,  assuming  an  identity  overlap  matrix. 


det|H"-^^-el|  = 0. 


(2.34) 


The  matrix  elements  or  integrals  are  the  key  assumptions  of  the  HMO,  where 


= {MH”I\<I,,)  = a. 


(2.35) 
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= (3,  and  u bonded 


0,  jjL  and  u not  bonded. 


(2.36) 


Thus  the  interelectronic  repulsions  are  averaged  into  unspecified  one-electron  integrals 
’Coulomb’,  a,  and  ’resonance’,  and  the  construction  of  the  Hamiltonian  matrix  is 
based  on  molecular  topology. 

Later  the  extended  Hiickel  theory  (EHT)  [14]  was  developed.  EHT  is  an  all- 
valence one-electron  theory.  The  idea  was  first  put  forth  by  Mulliken  in  1940s,  then 
further  developed  by  Wolfsberg  and  Helmholtz,  Lohr  and  Lipscomb,  and  Hoffman. 
The  one-electron  Hamiltonian  is  defined  as  the  following: 


where  is  the  valence  shell  ionization  potential  of  atomic  orbital  (j)^,  K is  a numerical 
constant  with  values  between  1.75  and  3 in  applications.  Unlike  the  HMO,  EHT  is 
a method  that  includes  overlap,  so  one  needs  to  solve  the  secular  equation  with  a 
non-identity  overlap  matrix. 


(2.38) 


(2.37) 


det(H"-^^  - Asi)  = 0. 


(2.39) 


Similar  to  the  HMO  method,  the  early  version  of  EHT  is  not  an  interative  procedure 
either  (however,  iterative  procedures  have  been  introduced  to  extend  the  usefulness 
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of  this  method).  These  two  methods  are  rather  valuable  for  conjugated  organic 
compounds.  However,  for  nonalternant  hydrocarbons,  or  many  systems  containing 
heteroatoms,  these  methods  do  not  even  offer  a qualitative  description. 

2.2.2  Semiempirical  SCF  MO  Methods  of  the  ZDQ  Type 

A 

Before  going  into  details  of  the  ZDO  type  methods,  we  would  like  to  have  our 
notations  defined  as  the  following:  F for  Fock  matrix,  and  F = H + G where  H is 
the  one-electron  matrix  and  G is  the  two-electron  matrix.  G =:=  J - |K  (close-shell), 
where  J is  the  Coulomb  integral  matrix  and  K is  the  exchange  integral  matrix. 

Based  on  the  resolution  of  the  HF  equations,  Pariser,  Parr  and  Pople  devel- 
oped a more  elaborate  method  (PPP)  [15-18].  This  method  is  an  SCF  method  but 
still  belongs  to  the  7r-electron  type.  Interelectronic  repulsions  are  included  in  the 
Hamiltonian  as  in  the  Hartree-Fock  method,  yet  similar  to  HMO,  overlap  is  ignored, 
and  further,  the  zero  differential  overlap  (ZDO)  approximation  is  made  for  the  two 
electron  integrals: 

= 0,  for  (2.40) 

which  results  in 

{^lX\v(7)  = (2.41) 

where  7^;^  are  often  treated  as  empirical  parameters.  Although  the  PPP  method 
is  not  really  applicable  to  systems  other  than  conjugated  and  aromatic  systems, 
it  is  a milestone  in  the  development  of  quantum  chemistry,  since  the  inclusion  of 
the  interelectronic  interactions  is  the  first  consistent  improvement  over  the  HMO 
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method.  Furthermore  the  ZDO  approximation  is  one  of  the  most  important  features 
for  subsequent  semiempirical  SCF  MO  methods. 

During  the  1960s,  Pople  and  collaborators  introduced  a series  of  ZDO  methods: 
the  complete  neglect  of  differential  overlap  (CNDO)  method,  intermediate  neglect  of 
differential  overlap  (INDO)  method  [19-22],  and  the  neglect  of  diatomic  differential 
overlap  (NDDO)  method.  These  are  all-valence-electron  methods,  with  the  overlap 
ignored  as  in  PPP  method.  However,  core  integrals  of  the  type 

(2.42) 

between  non-bonded  atoms  are  no  longer  ignored,  as  many  are  in  the  HMO  method. 
They  may  be  treated  in  a semiempirical  manner  in  order  to  accomodate  the  possible 
bonding  effect  of  overlap.  The  various  ZDO  methods  outlined  in  the  following  text  dif- 
fer mainly  in  the  extent  to  which  the  ZDO  approximation  is  invoked  in  the  calculation 
of  inter-electronic  repulsion  integrals.  Although  the  particular  ZDO  approximation 
used  will  be  different  among  the  various  ZDO  methods,  the  SCF  procedure  is  the 
same  for  all:  that  is,  the  orbital  energies  and  corresponding  molecular  orbitals  are 
obtained  by  solving  the  Roothaan  equation  self  consistently,  with  the  overlap  set  to 
the  identity  matrix. 

Before  discussing  the  various  methods  based  on  the  ZDO  approximation,  there 
is  a condition  that  all  these  methods  need  to  fulfill.  They  all  need  to  be  rotationally 
invariant  in  order  to  avoid  unphysical  results,  e.g.,  the  total  energy  varies  as  the 
coordinate  system  changes.  This  puts  certain  restriction  on  the  matrix  elements. 
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Complete  Neglect  of  Differential  Overlap  method 

As  indicated  in  the  name,  complete  neglect  of  differential  overlap  (CNDO),  the 
differential  overlap  considered  in  this  method  follows  strictly  Eqn.  (2.40).  From 
this  equation  it  follows  that  the  two-electron  integrals  in  the  CNDO  scheme  are 
approximated  as 

(/rAlr^Cr)  1 AB^nv^Xai  (2.43) 

where  orbitals  and  (j)y  are  centered  on  atom  A and  (j)x,  (j)a  are  centered  on  atom 
B.  Then  the  Fock  matrix  elements  are  given  as  the  following: 


F^lA^lA  — ~ ’^P^.a/j-aIaa  + '^PebIab,  (2.44) 

^ B 

PfJ-Al^A  ~l^PtXA^'A'yAA-i  (2.45) 

PfiA<^B  ~ PflA(TB  ~ l^PflAiyBlABl  (2-46) 

Pbb  = Y.P<r<r-  (2.47) 

<t€B 


In  these  equations  7^5  are  the  Coulomb  type  integrals,  which  depends  on  the  particu- 
lar pair  of  atoms  that  are  involved.  745  represents  the  average  electrostatic  repulsion 
between  any  electron  on  atom  A and  any  electron  on  atom  B,  and  it  is  made  to 
be  approximately  equal  to  as  the  interatomic  distance  Rab  increases.  Making 
these  two-center  two-electron  integrals  depend  on  the  nature  of  atoms  A and  B,  and 
Rab,  maintains  rotational  invariance.  The  core  operator  H is  defined 


(2.48) 
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where  -Vb  is  the  potential  for  an  electron  on  atom  A,  due  to  the  attraction  from 
nucleus  B,  and  the  repulsion  from  the  inner  shell  electrons  of  atom  B.  The  matrix 
elements  of  Eqn.(2.48)  in  the  given  atomic  orbitals  are: 


H, 


fJ-AtJ-A 


(/^a|  - -V^  - VaIha)  - ^ (ij.a\Vb\ha) 

UiiAtlA  ~ 

B^A 


^^^A‘'A 


— — 0, 

B^A 

= {^J■A\  — -V^  — Va  — VbIctb)  — ^ 

^ C^A,B 

— f^flAf^B  ~ PaB^I>'A<Jb- 


(2.49) 


(2.50) 

(2.51) 

(2.52) 


The  text  below  explains  the  approximations  for  Eqs.  (2.49)  - (2.52).  The 
term  in  Eqn.  (2.49)  is  a one-center  term,  and  it  contains  the  one-center  kinetic  en- 
ergy and  potential  energy.  Here  the  term  ’potential  energy’  refers  to  the  attraction 
with  the  nucleus  A on  which  (j)^  is  centered  and  the  repulsion  with  the  core  elec- 
trons centered  on  A.  This  term  vanishes  in  Eqn.  (2.50)  for  the  following  reason: 
within  the  given  atomic  orbitals,  the  one-center  term  becomes  zero  because 

the  valence  atomic  orbitals  and  (jy^,  belonging  to  the  same  atom  are  orthogonal 


to  each  other.  The  two-center  term  in  Eqn.  (2.50),  (/i^lVeli^a),  is  also  zero,  as  the 
monoatomic  differential  overlap  is  zero  according  to  Eqn.  (2.40).  Finally, 

the  term  Y.ci^A,B{AA\Vc\o' b)  in  the  two-center  one-electron  elements  in  Eqn.  (2.51)  is 
considered  comparable  to  the  three-center  two-electron  terms  and  should  be  omitted. 
However,  the  element  {ha\  ~ — Va  — Vb\(J b)  which  connects  two  atomic  centers  A 
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and  B and  therefore  governs  the  bonding  between  them,  is  most  commonly  approx- 
imated with  the  parameter  the  so  called  ’resonance  integral’.  /?  is  chosen  as 

an  atomic  parameter  in  order  to  satisfy  invariance  conditions,  and  there  are  different 
choices  of  the  forms  of  One  common  choice  is  given  by  the  Wolfsberg-Helmholtz 
formula  [23]: 

= + (2.53) 

Combining  the  core  integrals  in  Eqs.  (2.49)  - (2.52)  with  the  Fock  matrix  ele- 
ments, Eqs.  (2.44)  - (2.46)  gives 


F 

■*  tiAfJ-A 

— UfiAHA  + {Paa  ~ -zP ua^a)'7aa  + X (PbbJab  — Vab), 

^ B^A 

— ~'^PiiAi'A'y  aa, 

(2.54) 

F 

(2.55) 

F 

^l■A<^B 

~ !^AB^ilA<^B  ~ 2^t^A<^B'yAB-, 

(2.56) 

where 

Vab  = {^iA\VB\^^A)■  (2.57) 

Let  the  net  charge  on  atom  B be  Qb,  then 


Qb  = Zb-  Pbb,  (2.58) 


and  applying  this  charge  to  expression  of  gives 

— PfiAiJ-A  + [Paa  — -zPtiAHA)lAA  + X]  [~Qb1ab  + {ZbIab  ~ Vab)]-  (2.59) 

Bi^A 
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A widely  used  approximation  which  can  be  understood  from  the  rearranged  expres- 
sion above,  was  introduced  to  the  two-center  Fock  matrix  element.  The  first  part, 
Qb"7ab,  respresents  the  potential  that  arises  from  the  total  charge  on  atom  B,  and 
the  second  part, 

ZbIab  — Vab,  (2.60) 

is  referred  as  the  penetration  integral  following  Goeppert-Mayer  and  Sklar  [24].  In 
most  of  the  parameterization  schemes,  the  penetration  integrals  are  normally  consid- 
ered to  be  zero  in  order  to  compensate  for  corrections  that  are  needed  in  the  ZDO 
approximation,  and  hence  this  leads  to 

— ^tiAfJ-A  + {PaA  — -P UAiiA)lAA  + X]  {PbB  — ZB)'yAB-  (2-61) 

There  are  many  different  parameterization  schemes  used  with  the  CNDO  Hamil- 
tonian, but  these  will  not  be  reviewed  here.  However,  due  to  the  similarity  between 
CNDO  and  INDO,  some  of  the  details  will  be  mentioned  later  together  with  those 
schemes  for  the  INDO  method. 

Intermediate  Neglect  of  Differential  Overlap  method 

The  Hamiltonian  of  the  intermediate  neglect  of  differential  overlap  (INDO) 
method  is  similar  to  the  CNDO  Hamiltonian.  The  major  difference  is  that  one-center 
two-electron  integrals  of  the  type  {iiA<yAWA>^A)  are  retained. 
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The  one-center  one-electron  integral  elements  of  the  INDO  Hamiltonian  are 
~ (a*a|  — ^ {^J‘A\yB\^^A) 

^ B^A 

= - L ^ab,  (2.62) 

B^A 

where  (j)^  is  centered  on  atom  A.  is  the  atomic  core  integral,  which  stands  for  the 
one-center  kinetic  energy,  attraction  with  the  nucleus,  and  repulsion  with  the  core 
electrons  on  the  same  center,  in  the  same  way  as  it  is  defined  in  the  CNDO  method, 
similarly  for  the  term  Vab-  Also  the  two-center  core  integrals  are  approximated  the 
same  way,  and  if  the  Wolfsberg-Helmholtz  formula  [23]  is  adopted  one  obtains: 

^^^Al'B  ~ + PB)S^^,yg,  (2.63) 

where  Pa,  Pb  are  empirical  parameters  for  atom  A and  atom  B respectively. 

The  diagonal  one-center  two-electron  integral  elements  are 

Gfij^iiA  — X]  P\A<rA\{^^A^A\^J'A<^ a)  ~ 'A{^J‘A^A\(^ A^J‘A)\  + ^ PbbJAB,  (2-64) 

B^A 

where 

Pbb  = X! 

X€B 


(2.65) 


27 


It  is  important  that  the  parameterization  of  74^  satisfies  the  rotational  invariance 
condition.  In  a basis  set  with  s-  and  p-  functions  only,  Eqn.  (2.64)  becomes 

E P>^A>^A[{^^A^A\|J-A^A)  — -^{^^A^A\^A^J‘A)\  + X]  ^BBIAB-  (2.66) 

>^A  ^ B+A 

The  off-diagonal  one-center  two-electron  integral  elements  are, 

^tiAl'A  ~ E P\a<ta[{^^A^a\^A(^a)  — '^{^A^a\(^aI'a)\-, 

\?>{PaI'a\i^aIJ-a)  - {^J■Al^A\|J■AyA)\,  (2.67) 

where  the  second  equality  holds  when  only  s-  and  p-  types  orbitals  are  involved. 
Many  of  these  integrals  are  absent  in  CNDO  and  thus  constitutes  the  main  difference 
between  CNDO  and  INDO. 

Finally,  there  are  the  two-center  two-electron  integral  elements,  which  are  ex- 
change type  of  integrals  because  of  the  ZDO  approximation; 

~ ~'^Pb'A(^b^abi  (2.68) 

where  orbital  is  centered  on  atom  A and  is  centered  on  atom  B.  Combining 
both  the  one-electron  and  the  two-electron  matrix  elements  we  get  the  Fock  matrix 
elements  for  the  INDO  method  (with  s-  and  p-  orbitals  only  in  the  basis  set): 

^BA/^A  = ^BABA  + '^P\A>^A[{^^A^A\^^A>^A)  - -z{PA^A\^Afi-A)\ 

A4  ^ 
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+ X {PbbJab  — Vab), 

(2.69) 

B^A 

FfXAl^A 

= 2^tiA<'A[Hf^At^AWAf^A)  - if^Al^Alf^Al^A)], 

(2.70) 

F 

^AB^fiAO^B  ~ ’^FflAiyB'^ AB- 

(2.71) 

Again,  following  the  Goeppert-Mayer-Sklar  approximation  by  using  Eqn.  (2.58),  one 
can  rewrite  the  two  center  part  of  into  the  penetration  integrals,  Zb  Jab  ~ Vab, 
and  —QbJab,  and  ignore  the  penetration  integrals.  This  is  done  in  most  of  the  INDO 
schemes  and  leads  to  an  expression  similar  to  Eqn.  (2.61): 

F^aha  ~ UtJ.At>-A  PxaXa  \(f^A^A\l^A^A)  ~ ^((J-A^aI^AUa)] 

+ XI  (Pbb  - Zb)jab,  (2.72) 

B^A 

when  using  the  basis  set  consisting  of  s-  and  p-  functions  only. 

Neglect  of  Diatomic  Differential  Overlap  method 

In  the  NDDO  method  all  one-center  differential  overlaps  are  included;  viz., 

= ^4B0^^(1)0<t4(1)c^Wi-  (2.73) 

With  the  understanding  that  the  main  difference  is  in  the  two-electron  integrals,  we 
can  write  the  NDDO  Eock  matrix  elements: 

^fJ^AfJ^A  ~ ^t^AilA  ~ X {^^A\VB\pA)  + X Pl'AOA[{^^A^A\|J'A<yA)  (2.74) 

B^A  VAOA 
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~ c){^^A^AWAP'A)]-^  XI  PBB{^^A(^B\^J'A^B)1  (2.75) 

^ Bi^A 

= - X + X P‘^A>^Ai{f^A(7AWA>^A)  - 7:if^A(^A\>^Al^A)) 

B^A  (ta>^a  ^ 

“t”  ^ !!  F(jB^B^t^A(^ b\^a^b^  1 for  ^ 4^t^A  (2.76) 

B^A 

Fha(tb  ~ pAB^fiACTB  ~ 9 X Fi'A>^B(f^A(^BWA^B)-  (2.77) 

t'A^B 

The  potential  (^^Ihfikyi)  is  divided  into  four  groups  and  calculated  according  to  the 
Goeppert-Mayer-Sklar  approximation  [24]; 


(s/i|Vb|s4) 

~ Z*b{^  AS  b\s  aS  b)  1 

(2.78) 

{saIVbIpcta) 

= Z*ji{sASB\p<yASB), 

(2.79) 

{p(yA\VB\p(yA) 

= Z*B{paASB\p(TASB), 

(2.80) 

{pt^a\Vb\pt^a) 

= Zg(j}'KASB\pT^ASB)- 

(2.81) 

Z*g  in  the  expressions  above  is  the  effective  core  charge,  which  is  obtained  by  sub- 
tracting the  number  of  the  core  electrons  from  the  atomic  number  in  order  to  be 
consistent  with  the  valence-electrons  only  method,  and  reflects  the  effect  of  screening 
the  nucleus  from  the  valence  electrons  by  the  core  ones. 

The  two-center  two-electron  integrals  in  Eqn.  (2.77)  are  particular  for  the 
NDDO  method,  {ij-a^bWa^b),  and  they  can  be  evaluated  theoretically  using  Slater 
type  orbitals  or  by  a parameterization  scheme  by  Dewar,  the  modified-NDDO  (MNDO) 
method  [25, 26].  In  this  case  the  integrals  are  obtained  by  doing  a multipole-multipole 
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expansion  of  the  corresponding  integrals  [27,  28]: 

{^iA(rBWAXB)  = E (2.82) 

h,l2,m 

where  M;  is  the  2^-pole  moment  of  the  charge  distribution  arising  from  the  corre- 
sponding atomic  orbital  product,  and  the  subscripts  /,  and  m correspond  to  the  order 
and  the  orientation  of  the  given  multipole.  Each  multipole  can  be  represented  by  an 
appropriate  configuration  of  2*  point  charges,  the  multipole  interaction  energy  can 
be  calculated  by  the  following  function, 

1 2l\  2/2 

^ E 04-  + [(a1  + (2.83) 

^ i=l j=l 

where  the  parameters  and  are  set  by  the  limit  at  Rab  = 0.  Clearly  as  Rab  ^ 
00,  the  function  approaches  1/Rab 

2.2.3  Parameterization  for  the  INDQ  Method 

The  main  concern  in  this  work  is  to  find  a method  useful  for  both  geometry  and 
electronic  excitation  calculations.  The  general  INDO  Hamiltonian  has  been  outlined 
earlier  in  the  text.  However,  there  are  many  versions  of  the  INDO  method  all  using 
different  parameterization  schemes.  It  is  appropriate  to  outline  briefly  the  most 
important  parameterization  schemes  before  I proceed  to  discuss  my  own  work. 
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INDQ  parameterization 

The  INDO  parameterization  is  similar  to  the  scheme  for  CNDO/1  or  CNDO/2. 
The  integrals  that  are  parameterized  are  the  overlap  integrals,  5^^,  the  core  Hamilto- 
nian integrals  Vabi  the  resonance  integrals  and  the  two-electron  repulsion 
integrals  74^  and  7^^.  For  the  INDO  model,  also  all  the  one-center  two-electron 
integrals  are  required,  which  is  the  difference  compared  to  CNDO  model.  The  basis 
set  for  the  valence  shell  consists  of  Slater  type  orbitals  [29]: 

where  ^ is  the  orbital  exponent  calculated  according  to  Slater  rules  [29]  (except  for 
hydrogen,  where  ^=1.2  is  used  as  this  value  is  more  appropriate  for  hydrogen  in 
molecules),  Oq  is  the  Bohr  radius,  n is  the  principle  quantum  number,  I is  the  angular 
quantum  number,  m is  the  magnetic  quantum  number,  Yi^{6,(f))  are  the  spherical 
harmonics,  and  {r,9,<j))  are  the  spherical  polar  coordinates.  The  one-center  basis 
functions  are  orthogonal  to  each  other,  except,  for  example,  the  Slater  s functions 
from  different  shells  are  not  orthogonal  to  each  other  because  of  the  lack  of  the  nodes 
in  the  radial  part.  However,  since  only  the  valence  shell  is  considered  in  the  all- 
valence semiempirical  methods,  no  problem  arises  from  this  feature.  The  overlap 
integrals  are  explicitely  calculated.  The  two-center  two-electron  integral  7^b,  which 
is  treated  as  an  average  interaction  between  electrons  in  the  valence  orbitals  on  atom 
A and  B,  is  evaluated  as  two-center  Coulomb  integrals  using  the  s-orbital  but  use 
the  exponents  of  the  corresponding  orbitals  (as  was  done  in  SINDOl  [30])  in  order 
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to  preserve  the  radial  extent: 


(2.85) 


and  the  one-center  two-electron  integrals  are  also  calculated  ab  initio  too. 

The  potential  energy  involving  a valence  electron  on  atom  A with  the  core 
electrons  and  the  nucleus  of  atom  B,  is  calculated  from  the  following  equation 
in  the  CNDO/1  scheme, 


where  Zb  is  the  core  charge  of  atom  B and  Tib  is  the  distance  of  electron  1 from  the 
nucleus  B.  An  effective  core  charge  Z*g  is  introduced: 


where  is  the  number  of  core  electrons  on  atom  B.  The  effective  core  charge 
provides  the  screening  effect  exerted  by  the  core  electrons  so  that  the  valence  elec- 
trons feel  the  attraction  from  a nucleus  with  only  the  effective  core  charge.  We 
replace  the  core  charge  in  the  calculation  of  Vab  with  the  effective  core  charge.  The 
equilibrium  bond  lengths  obtained  from  the  CNDO/1  method  tend  to  be  too  short, 
and  the  dissociation  energies  tend  to  be  too  large.  The  error  comes  partly  from  the 
penetration  integrals,  and  partly  from  the  nature  of  the  ZDO  approximation.  The 


f 4(1)— 

J riR 


(2.86) 


Zg  = Zb  ~ n\ 


core 


(2.87) 
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problem  can  be  resolved  by  ignoring  the  penetration  integrals,  i.e.,  according  to  the 
Goeppert-Mayer-Sklar  approximation  [24],  the  evaluation  of  becomes 

Vab  - {i^a\Vb\^ia)  = Z*b^ab-  (2.88) 

This  is  implemented  in  the  CNDO/2,  INDO/1  and  INDO/S  parameterizations  as 
well. 

For  the  resonance  integral  /?%,  the  Wolfsberg-Helmholtz  formula  [23]  is  applied, 

(2.89) 

where  /?/  is  an  empirical  parameter  specific  for  atom  I. 

The  core  integrals  contain  the  one-center  kinetic  energy  and  the  core  po- 
tential. These  integrals  are  obtained  from  the  observed  atomic  energy  levels.  For  an 
atom  A with  electron  configuration  (2s)"* (2p)”,  the  total  energy  within  the  CNDO 
approximation  is  given  by  [19] 

E(X,  2s"*2p")  = mU2s2s  + nll2p2p  + — "^—^^yAA,  (2.90) 

2!(m-|-n  — 2)! 

where  X is  the  filled  core.  When  one  p-electron  is  taken  away,  the  total  energy 
becomes 

£(.Y,2»’”2p”-‘)  = mU,.,.  + (n  - 

2!(m  -I-  n — 3)! 


(2.91) 
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Ignoring  orbital  relaxation  after  the  loss  of  one  electron,  the  ionization  potential 
derived  from  this  scheme  is 


and  therefore  the  core  integrals  U2s2s  and  C/2p2p  are  obtained  by  the  expression  below: 


Using  merely  the  ionization  energy  to  calculate  the  core  integral  is  the 
strategy  adopted  for  the  CNDO/1,  INDO/2,  and  INDO/S  methods.  However,  in 
order  to  account  for  the  ability  of  an  atomic  orbital  to  both  aquire  and  lose  electrons, 
both  the  ionization  potential  and  the  electron  affinity  are  used  in  CNDO/2  and  are 
related  to  through 


Ij,  = E{X,2s^2p^-^)  - E{X,2s^2p^) 


U2p2p  -{m  + n-  1)7.44, 


(2.92) 


— ^flA  i^A  ~ 1)744- 


(2.93) 


(2.94) 


which  is  obtained  by  combining  Eqn.  (2.93)  and  the  following: 


(2.95) 


where  is  the  electron  affinity  of  the  atomic  orbital  on  atom  A.  This  way 
of  employing  the  ionization  potential  and  electron  affinity  into  the  determination  of 


35 


core  integrals  assumes  frozen  orbitals  upon  losing  or  obtaining  an  electron,  as  shown 
in  Eqn.  (2.94).  As  discussed  previously  Koopmans’  theorem  is  a good  approximation 
for  ionization  potential  but  not  for  the  electron  affinity.  Therefore,  it  may  not  be  a 
surprise  that  CNDO/2  does  not  predict  ionization  potentials  very  well.  For  the  one- 
center  two-electron  integrals  in  the  INDO  method,  the  Slater-Condon  parameters 
[31,  32]  and  which  are  two-electron  integrals  involving  the  radial  parts  of 

the  atomic  orbitals,  are  used  to  evaluate  the  related  integrals: 


(ss|ss) 

= (spxlspx)  = F^  = 7aa, 

(2.96) 

(ssIpxPx) 

- 

(2.97) 

(PxPxlPyPy) 

II 

(2.98) 

{PxPx  \PxPx) 

II 

O 

+ 

(2.99) 

(PxPylPxPy) 

1 

II 

(2.100) 

* 


These  one-center  two-electron  terms  in  the  INDO  method  make  the  evaluation  of 
core  integrals  more  complicated  than  in  CNDO,  since  in  the  INDO  scheme,  the 
total  energy  for  an  atom  A with  electron  configuration  (2s)”"  (2p)"  becomes  [21] 

£(,Y,  2»”>2p“)  = mU,„.  + nU,,,,  + <!!L+  "X”  + ” ~ 

- ~mnG^  - -^n(n  - 1)F^. 

6 25 


(2.101) 
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For  example,  if  similar  parameterization  scheme  as  CNDO/2  is  followed,  then  for  the 
elements  B to  F in  the  second  row,  U2p,2p  is  given  within  INDO  (i.e.,  INDO/2)  by 


U, 


2p2p 


{hp  + ^2p) 


- (Z* 


mn  + lG\  + ^^{Z--5/2]Fl 


2 


(2.102) 


where  Z*  — m + n is  the  effective  core  charge. 

Within  the  Born-Oppenheimer  approximation,  the  nuclear-nuclear  repulsion  po- 
tential needs  to  be  added  to  the  electronic  energy  to  obtain  the  total  energy.  For 
the  all-valence  methods,  the  compatible  nuclear-nuclear  repulsion  is  evaluated  using 
the  effective  core  charge.  This  can  be  illustrated  using  the  case  of  a simple  diatomic 
molecule  with  the  numbers  of  core  electrons  being  and  ns-  We  assume  that  the 
core  electrons  are  rather  localized  and  hence  the  electron  exchange  between  the  cores 
are  ignored: 


E„ 


Z>aZb  _ ^aZb 
Rab  Rab 

Rab 


^bZa  riAriB 
Rab  Rab 


(2.103) 


The  total  energy  is  thus  calculated  from 


F, 


7*  7* 

— T?  I ^A^B 

total  ^elec  ' / . jy 

A>B 


(2.104) 
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INDO/S  and  CNDQ/S  parameterization 

Inspired  by  an  observation  of  Pariser’s  [18]  concerning  the  two-electron  integrals 
in  the  PPP  method,  Jaffe  and  Del  Bene  modified  the  CNDO/2  parameters  [33^35], 
which  made  it  possible  to  utilize  the  semiempirical  method  for  electronic  spectra 
calculations.  According  to  Pariser,  the  one-center  two-electron  integrals  are  given  by 


(2.105) 


where  is  the  electron  affinity  of  atomic  orbital  (f>^.  For  the  two-center  two-electron 
integrals,  at  first  they  used  the  Pariser-Parr  formula  [17]: 

7AS  = ^(74  + 7b)  - qRab  - bRlg,  (2.106) 


in  which  a,  and  b are  constants  that  are  determined  for  Rab  between  2.8  A and  3.7 
A.  Later  suggested  by  Zerner  and  Ridley  in  their  INDO/S  method  (see  below),  these 
integrals  were  evaluated  using  the  Mataga-Nishimoto  formula  [36]: 


Jab 


1 

a + Rab  ’ 


(2.107) 


where  a is  the  inverse  of  the  average  of  the  one-center  two-electron  integrals  of  atoms 
A,  and  B: 


2 

JHA  J<TB 


a = 


(2.108) 
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The  expression  of  the  resonance  integral  is  also  modified 

= + (2.109) 

where  the  overlap  matrix  element  is  modified  according  to 

I 

^0-/ct(/^(T4  |(7'o-b)  + (^TT  + 57r')/7r(/^7r^  Ic’^TTg)  + ■ • • , (2.110) 

where  the  factors  gi  are  the  Eulerian  rotational  factors  necessary  in  evaluating  the 
integrals  in  the  local  coordinate  system  and  transforming  back  to  the  molecular  co- 
ordinate system.  The  parameter  / is  to  remedy  the  problem  in  CNDO/2  of  overesti- 
mating the  interactions  among  tt  bonding  orbitals,  compared  to  a bonding  orbitals. 
The  optimized  values  for  / are; 

fa  = 1-0,  (2.111) 

U = 0.585.  (2.112) 

In  the  1970’s,  Zerner  and  Ridley  [37, 38]  developed  the  spectroscopy  version  of  INDO 
method,  called  INDO/S,  by  modifying  the  INDO  parameters.  Similar  to  what  was 
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done  in  CNDO/S,  the  one-center  two-electron  integrals  are  calculated  from  the  ex- 
perimental ionization  potentials  and  electron  affinities.  Different  weighting  factors 

Ms)  = 1.0,  (2.113) 

A = 0.585,  (2.114) 

Mp)  = 1.267,  (2.115) 


are  applied  to  the  overlap  between  different  types  of  molecular  orbitals  when  con- 
structing two-center  one-electron  integral  as  in  Eqn.  (2.109). 

The  Mataga-Nishimoto  form  of  the  two-center  two-electron  integrals  is  adopted, 
as  in  CNDO/S,  but  the  Weiss  parameter  fy  = 1.2  is  found  to  yield  better  results  for 
calculations  of  spectroscopic  properties  of  large  molecules: 

Moreover,  the  Slater-Condon  factors  are  from  atomic  spectroscopic  data  in  order  to 
be  compatible  with  the  one-center  7^^  integrals.  It  turns  out  that  the  experimental 
Slater-Condon  factors  are  about  0.6  times  the  corresponding  ab  initio  values.  The 
parameterization  of  INDO/S  was  carried  out  at  the  single-excitation  Cl  (CIS)  level 
and  based  on  experimental  spectra.  The  corresponding  nuclear-nuclear  repulsion 
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potential  for  calculating  the  total  energy  is  evaluated  in  terms  of  the  two-center  two- 
electron  integrals  and  the  effective  core  charges: 

Efiu—nu  — ^a^b^ab-  (2.117) 

A>B 

So  far  as  one  can  find  from  the  literature,  all  the  integrals  which  are  used  in  the 
spectroscopic  version  of  ZDO  methods  are  smaller  in  magnitude  than  those  which 
are  used  in  the  corresponding  ab  initio  case.  This  is  due  to  two  factors:  dynamic 
electron  correlation  and  the  effect  which  is  caused  by  using  the  minimal  basis  set 
in  the  model.  Including  electron  correlation  or  employing  big  basis  set  reduces  the 
effective  magnitude  of  the  ab  initio  integrals. 

A comparison  between  the  repulsion  using  the  Ohno-Klopman  formula,  the 
Coulomb  form  and  the  STO  two-electron  integral  (the  exponent  of  C is  set  to  1.625) 
is  shown  in  Fig.  2.1.  The  asymptotic  values  at  i?  = 0,  7,  in  the  case  of  the  C2 
molecule  is  indicated  in  the  figure. 
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Figure  2.1.  Comparison  between  Ohno-Klopman,  Coulomb  repulsion,  and  Slater 
two-electron  integral 


2.2.4  General  Remarks 

The  INDO  model  has  several  advantages  over  CNDO  in  that  the  INDO  method 
can  separate  n — > tt*  triplet  and  singlet  states  successfully  because  of  the  inclusion 
of  one-center  two-electron  integrals.  For  example,  the  singlet-singlet  excitation  from 
an  occupied  orbital  i to  a virtual  orbital  / is 


^Eif  — e f — £i  — Jif  -h  2Kif, 


(2.118) 
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and  the  excitation  energy  from  orbital  i to  / resulting  in  a triplet  state  is 

AE^f=ef-£i-Jif,  (2.119) 

where  J and  K are  the  molecular  Coulomb  and  exchange  integrals  respectively.  In  a 
CNDO/S  calculation,  these  two  types  of  excitations  are  degenerate  as  Kif  = Knn  = 0. 
In  addition,  CNDO  is  useless  in  the  calculation  of  transition  metal  systems  as  these 
neglected  exchange  integrals  are  often  of  the  order  of  an  atomic  unit  [39].  However, 
for  calculating  bond  lengths,  bond  angles,  and  electron  density  distributions,  INDO 
does  not  offer  too  much  of  an  improvment  over  CNDO.  NDDO,  on  the  other  hand, 
by  including  two-center  two-electron  integrals,  (/x^AbIz/^ctb),  has  the  potential  to  be 
superior  to  the  original  INDO  method  for  calculating  geometries  and  directional  prop- 
erties, and  thus  predicting,  e.g.,  molecular  dipole  moments.  It  is  not  clear,  however, 
whether  this  potential  is  realized  as  SINDO  [30],  an  INDO-based  model,  seems  as 
effective  as  AMI  [40]  (an  NDDO  method)  in  predicting  geometries  and  thermochem- 
istry. Although  many  different  parameterization  schemes  have  been  proposed  in  the 
past,  none  of  them  can  be  considered  accurate  enough  for  all  molecular  properties. 
That  is  to  say  that  either  the  choices  of  parameters  are  not  optimal,  or  the  model 
Hamiltonian  is  not  accurate,  and  therefore  the  applicability  of  current  semiempirical 
methods  often  depends  on  the  purpose  and  the  means  for  which  the  model  is  applied. 


43 


Table  2.1.  The  comparison  of  all- valence  electron  methods 


Method 

Geometry 

Charge  distribution 

Electronic  excitation 

CNDO/2 

good 

good 

poor 

CNDO/S 

not  app. 

poor 

good 

INDO 

good 

good 

poor 

INDO/S 

not  app. 

good 

good 

MINDO/3 

good 

good 

poor 

MNDO 

good 

good 

poor 

AMI  [40] 

good 

poor 

poor 

PM3  [41] 

good 

poor 

poor 

NOTE:  All  comments  are  from  [42],  except  for  CNDO/S,  INDO/S,  AMI,  and  PM3. 

It  is  clear  from  the  comments  on  the  performance  of  major  semiempirical  meth- 
ods that  not  many  of  them  are  good  for  calculating  electronic  spectra.  Furthermore, 
the  only  two  methods  which  are  good  for  calculating  spectra  are  bad  for  geometry  op- 
timization. Hence,  this  gives  a strong  motivation  to  the  author  to  accomplish  a new 
semiempirical  model  which  is  multi-purposed,  especially  equipped  with  the  ability  of 
calculating  both  electronic  excitation  and  geometry,  and  thus  applicable  to  problems 
in  photochemistry. 


CHAPTER  3 
XINDO 

3.1  Features 

The  XINDO  method  is  designed  to  be  applicable  to  both  geometry  and  electronic 
excitation  calculations.  A common  set  of  parameters  will  be  sought.  As  stated  in  the 
theory  review  in  the  previous  chapter,  current  existing  methods  are  not  able  to  yield 
both  excitation  energies  and  geometries  equally  well  with  a single  set  of  parameters. 
Among  those  methods,  the  INDO/S  method  has  been  used  to  calculate  molecular 
spectra  quite  successfully,  but  this  is  not  the  case  for  the  other  similar  method,  the 
INDO/1.  On  the  other  hand,  the  INDO/1  method  gives  good  molecular  geometries 
but  the  INDO/S  method  does  not.  The  main  difference  between  these  two  methods  is 
the  choice  of  two-electron  integrals.  Therefore,  a promising  approach  is  to  understand 
the  reason  why  the  INDO/S  method  is  not  good  for  calculating  geometries.  Then 
one  can  take  advantage  of  the  INDO/S  model  which  has  proven  good  for  electronic 
excitation  calculations,  and  adjust  it  to  also  yield  good  molecular  geometries. 

We  use  the  INDO/S  two-electron  integrals  for  the  XINDO  method,  since  the 
INDO/S  method  has  been  successful  in  electronic  excitation  calculations,  and  the 
key  reason  is  the  adoption  of  empirical  parameters  for  the  two-electron  integrals. 
However,  the  development  so  far,  only  applies  to  basis  set  with  s-  and  p-  type  orbitals. 
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With  INDO/S  as  the  starting  point,  the  modifications  listed  below  are  made  in  order 
to  obtain  a new  method  which  can  be  used  to  calculate  geometry: 

1.  penetration  integrals  (expression  (2.60))  are  included  by  using  a separate  param- 
eter for  evaluating  the  nuclear-electron  attraction. 

2.  an  effective  core  potential  (ECP)  [43]  is  added  to  provide  a short  range  repulsion 
between  core  electrons  and  valence  electrons.  This  repulsion  has  much  sharper  dis- 
tance dependence  at  short  interatomic  distance  than  the  overlap.  ECP  can  be  used 
to  compensate  the  possible  drawback  arising  from  the  inclusion  of  the  penetration 
integrals. 

3.  off-diagonal  kinetic  energy  elements  are  included.  Again,  the  distance  dependence 
is  also  very  different  than  the  overlap. 

4.  two-center  one-  and  two-electron  integrals  are  approximated  using  the  Mulliken 
approximation  [44]  together  with  the  diagonal  one-  and  two-electron  integrals,  re- 
spectively. 

5.  in  the  spirit  of  the  EHT,  XINDO  will  include  overlap  and  solve  the  secular  equation 
with  non-identity  overlap. 

6.  the  STO  Coulomb  matrix  is  constructed  with  the  LTO  density  matrix.  This  allows 
us  to  use  the  ZDO  approximation  (INDO  or  NDDO)  to  avoid  calculating  complex  3- 
and  4-centered  integrals. 

7.  the  nuclear-electron  integrals  are  used  for  evaluating  the  nuclear-nuclear  repulsion 


energy. 
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Besides  yielding  both  geometry  and  spectra,  there  are  some  other  useful  features 
arising  from  the  way  we  design  XINDO:  the  inclusion  of 

1.  overlap  improves  the  description  of  the  energies  of  the  antibonding  molecular 
orbitals, 

2.  multi-center  integrals  improve  the  description  of  conformational  properties,  dipole 
moment,  etc.. 

In  the  following  section,  the  details  of  the  work  and  the  rationale  of  the  modifi- 
cations will  be  given. 

3.2  Details  of  Building  XINDO 

3.2.1  Comparison  with  Previous  Models 

It  is  natural  to  look  into  the  ZDO  approximation  and  the  Hiickel  type  methods 
once  more.  Doing  so,  one  may  understand  what  causes  the  constraint  of  having  to 
choose  between  spectroscopic  properties  or  geometry  in  the  ZDO-based  methods,  and 
one  may  also  find  useful  approximations  for  our  method  development. 

EHT 

The  EHT  [14]  method  is  a one-electron  theory,  i.e.,  without  two-electron  inte- 
grals in  the  Fock  matrix.  Therefore  EHT  is  not  ideal  for  calculating  the  electrostatic 
properties.  However,  since  the  overlap  matrix  is  not  set  to  identity  in  EHT,  it  is  able 
to  account  for  the  shapes  of  orbitals,  and  hence  benefits  the  calculations  of  molec- 
ular geometries.  The  use  of  the  overlap  matrix  also  allows  the  unsymmetric  energy 
splitting  of  the  bonding  and  antibonding  MO’s  as  described  below. 
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Consider  a two-orbital  homonuclear  EHT  Hamiltonian  matrix: 


H = 


H\i  Hi2 
Hn  Hu 


(3.1) 


From  this  Hamiltonian  and  the  secular  equation  Eqn.  (2.18),  one  obtains  the  orbital 
energies, 


e = 


H\\  ± H\2 

1 ± Si2 


(3.2) 


If  the  ZDO  approximation  is  applied  in  the  above  expression,  S12  = 0,  then  the 
molecular  orbitals  are  formed  with  a symmetric  split  from  the  constituent  atomic 
orbitals,  which  makes  the  antibonding  orbital  higher  in  energy  than  the  bonding 
orbital  by  2\Hi2\.  Consider  now  the  situation  where  the  overlap  between  the  diffuse 
atomic  orbitals  approaches  1.  In  such  a case  the  antibonding  orbital  is  elevated  in 
energy  far  greater  than  2 1 i/12 1,  and  the  bonding  MO  is  depressed  to  a much  less 
degree.  This  can  be  understood  from  Eqn.  (3.2),  as  the  denominator  approaches 
zero,  the  antibonding  MO  energy  increases  dramatically.  The  illustration  given  here, 
and  shown  in  Fig.  3.1,  indicates  that  the  overlap  elements  should  be  included  to 
improve  on  the  energy  gap  between  bonding  and  antibonding  orbitals.  The  inclusion 
of  the  overlap  also  allows  the  eventual  use  of  better  than  a minimal  basis  set,  as 
one-center  overlap  is  also  included.  Therefore,  one  obvious  way  to  improve  on  the 
INDO  method  is  to  express  the  XINDO  Hamiltonian  and  the  parameterization  in 
a non-orthogonalized  basis  set,  e.g.,  the  STOs,  and  then  solve  the  secular  equation 
with  a non-identity  overlap  matrix. 
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ZDO  methods 

CNDO/S  and  INDO/S  methods  are  among  the  very  few  semiempirical  methods 
which  can  calculate  electronic  spectroscopic  properties.  However,  neither  can  be  used 
to  calculate  molecular  geometries.  As  mentioned  earlier,  the  CNDO/S  method  can- 
not calculate  transition  metal  systems  correctly  but  INDO/S  can,  also  the  INDO/S 
can  distinguish  the  singlet  to  singlet  state  electronic  excitation  from  the  singlet  to 
triplet  state  excitation  in  the  cases  of  cr  — y tt*  and  tt  — y a*  excitations.  Therefore 
we  will  deal  with  the  INDO/S  model  rather  than  CNDO.  Another  method  at  the 
INDO  level,  the  INDO/1  method,  gives  satisfactory  molecular  geometries,  but  not 
electronic  excitations.  The  main  difference  between  the  INDO/S  and  INDO/1  meth- 
ods are  the  two-electron  integrals.  The  two-electron  integrals  that  are  used  in  the 
INDO/1  method  are  evaluated  ab  initio  over  STOs;  the  ones  that  are  used  in  the 
INDO/S  method  are  evaluated  according  to  Pariser  [18]  from  atomic  experimental 
data.  Obviously  such  integrals  contain  electron  correlation  and  compensate  for  the 
effects  of  the  use  of  a minimal  basis  set,  and  thus  make  the  calculations  of  electronic 
spectra  possible.  However,  the  reason  why  INDO/S  cannot  be  used  to  calculate 
geometry  is  still  unclear. 

Recall  the  commonly  used  Goeppert-Mayer-Sklar  approximation  [24]  in  all  ZDO 
methods.  In  that  approximation,  the  penetration  integral  ZbJab  ~ Vab  is  neglected. 
Vab  denotes  the  net  attraction  between  an  electron  centered  on  atom  A and  nucleus 
B.  ZbJab  denotes  the  repulsion  between  the  same  electron  on  atom  A and  the 
electrons  on  atom  B.  Applying  this  approximation  means  that  one  calculates  the 
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nuclear-electron  attraction  by  using  the  two-electron  integral  'Jab-  And  in  the  case 
of  INDO/S,  this  approximation  imposes  the  same  extent  of  electron  correlation  and 
minimal  basis  set  effect  on  both  the  two-electron  integral  and  the  nuclear-electron 
integral.  An  unnessary  constraint  may  occur  from  such  assumption. 

To  date,  there  is  no  report  of  electronic  excitation  calculations  using  the  NDDO 
method.  However,  this  method  may  have  some  advantages  in  geometry  calculation 
because  it  contains  the  two-center  two-electron  integrals  {^XAcyB\^A^B)-  The  only  two- 
center  two-electron  integrals  existing  in  the  INDO  methods  are  of  the  exchange  type. 
The  fact  that  the  INDO  methods  tend  to  overestimate  the  stabilization  of  ring  forms 
over  bent  or  linear  forms  for  certain  molecules  indicates  that  it  may  be  necessary  to 
include  these  two-center  two-electron  Coulomb  integrals  to  remedy  this  deficiency  in 
the  INDO  methods.  However,  we  note  the  success  of  the  SINDO  model  which  is  an 
INDO  model,  with  reported  equal  accuracy  to  the  AMI  and  PM3  methods  and  both 
are  NDDO  models. 

3.2.2  Integral  by  Integral 

In  this  section,  all  integrals  that  enter  the  XINDO  Fock  matrix  will  be  discussed. 
The  XINDO  method  is  constructed  in  the  STO  representation,  which  is  different  from 
all  other  ZDO  methods.  Therefore,  we  want  to  know  how  the  same  Hamiltonian 
differs  in  the  STO  representation  and  in  the  ZDO-compatible  representation. 

The  ZDO  approximation  is  often  justified  by  expressing  the  Hamiltonian  in 
the  Lowdin  symmetric  orthogonalized  basis  [45,46].  We  therefore  will  study  the 
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differences  with  the  one-electron  and  two-electron  integrals  in  these  two  different 
basis  sets,  STOs  and  LTOs. 

The  Lowdin  type  orbitals  have  the  important  property  that 

(Xf  (3.3) 

and  yet  these  orbitals  resemble  the  orginal  nonorthogonal  atomic  orbitals  as  much  as 
possible  in  the  least-square  sense,  i.e., 

- X^l)dv\  = min.  (3.4) 

One  may  transform  STOs  to  LTOs  by  the  following  expression: 

^ ^sro^-i/2^ 

where  A = 1 -f-  S is  the  overlap  matrix  of  the  STOs.  It  is  clear  that  the  differential 
integral  matrix  between  the  LTOs  is  the  identity  matrix. 

^LTOt^LTO  ^ ^~^/2^ST0-^^ST0^-1/2 

= 1,  (3.6) 

which  is  consistent  with  one  of  the  ZDO  approximations.  As  far  as  the  two-electron 
integrals  which  are  retained  in  the  ZDO  schemes  are  concerned,  no  drastic  change 
shows  up  when  they  undergo  symmertic  orthogonalization  transformation.  It  also 
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seems  a reasonable  assumption  that  by  applying  empirical  parameters  to  these  in- 
tegrals, the  small  changes  introduced  by  the  transformation  from  the  STO  to  LTO 
basis  may  be  taken  care  of.  Moreover,  with  a symmetrically  orthogonalized  basis  set, 
the  electron-repulsion  integrals  neglected  in  the  ZDO  schemes  are  generally  found  to 
be  quite  small  [47]. 

The  effects  caused  by  orthogonalization  transformation  of  the  one-electron  Hamil- 
tonian matrix  are  illustrated  by  the  following  expressions. 

^LTO  ^ 

which,  in  the  case  of  a two  orbitals  of,  for  example,  a diatomic  molecule,  leads  to  the 
following: 


ttLTO 

= - si.) 

(3.8) 

itLTO 

(3.9) 

= + HlJ°) 

(3.10) 

ttLTO 

_ zjSTO  q ttLTO 

(3.11) 

where  orbital  is  on  atom  A,  and  (j)^  is  on  atom  B.  Assume  that  the  difference  be- 
tween and  introduced  by  this  orbital  transformation  is  the  orthogonal- 

ization correction.  Numerical  results  show  that  such  correction  for  diagonal  elements 
1-^/^'™  ~ is  small,  but  the  off  diagonal  elements  and  are 

of  similar  magnitude,  since  <<  [48].  Thus  one  may  conclude  that  the 
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ofF-diagonal  one-electron  matrix  elements  of  the  XINDO  are  expected  to  be  larger 
than  the  ones  in  INDO/S,  but  the  diagonal  elements  should  be  similar  in  STO  or 
LTO. 

We  now  turn  our  attention  to  the  two-electron  integrals.  The  change  in  one-  and 
two-center  integrals  is  small  when  the  ZDO  two-electron  integrals  are  transformed 
from  the  STO  to  LTO,  whereas  integrals  involving  differential  diatomic  overlap  be- 
come quite  small.  Fischer- Hjalmars  [49],  and  more  recently,  Thiel  and  Weber  [50] 
have  shown  that  after  applying  an  expansion  of  the  overlap  matrix 

A-'/"  = (1  + = 1 - Is  + + ■ . . , (3.12) 

to  the  right  of  the  following  expression, 

f^LTO  ^ ^-1/2 STO ^-1/2^  ^3  ^3^ 

and  collecting  terms  up  to  first  order  in  S,  the  transformed  integrals  are  exactly  the 
same  as  in  the  non-orthogonal  basis  set  at  NDDO  level.  If  the  collected  terms  are  up 
to  second  order  in  S,  then  in  addition  to  the  NDDO  type  integrals  some  very  small 
two-center  integrals  are  introduced.  Based  on  this  observation,  no  obvious  change  in 
the  two-electron  integrals  is  expected.  Similarly,  one  will  not  expect  very  different 
two-electron  integrals  in  the  XINDO  method  than  other  ZDO  methods. 
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Nuclear-electron  integral 

From  the  conclusion  after  reviewing  the  ZDO  methods,  we  decide  to  abandon 
the  Goeppert-Mayer-Sklar  approximation  [24],  In  order  to  do  so,  we  will  introduce 
a parameter  for  the  one-center  nuclear-electron  attraction  = 7^^.  This 

'y^lA  calculate  the  two-center  nuclear-electron  attractions.  In  other 

words,  we  no  longer  borrow  the  two-center  two-electron  integrals  to  evaluate  the 
two-center  nuclear-electron  attraction.  Coffey  [51]  showed  that  ignoring  penetration 
integrals  was  equivalent  to  ignoring  an  effective  core  potential  and  orthogonalization 
corrections,  which  are  terms  that  we  are  willing  to  include  explicitly. 

The  use  of  different  weighting  factors  for  the  two-center  one-electron  integrals 
in  the  INDO/S  method,  Eqs.  (2.113)  - (2.115),  indicates  that  the  parameters  for  tt 
and  a bonds  should  be  different.  In  the  XINDO,  the  nuclear-electron  integral  7^^  is 
an  orbital  dependent  atomic  parameter. 

We  then  proceed  with  nuclear-electron  attraction  between  different  atoms  Vab- 
When  treating  two-center  terms  like  (/i/i|VB|//A)  = ^^e  most  important  con- 

cern is  to  have  a proper  distance  dependence  in  the  proposed  functional  form  of 
the  parameterization  that  can  represent  the  two-center  interactions  above.  In  the 
INDO/S  model  either  the  Mataga-Nishimoto  form  [36]  as  in  Eqn.  (2.107)  or  the 
Ohno-Klopman  form  [52,  53] 


7, 


(4(: 


+ 


(3.14) 
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is  used  for  the  two-electron  interactions,  and  as  a result,  good  spectra  can  be  obtained. 
The  Mataga-Nishimoto  (MN)  and  the  Ohno-Klopman  (OK)  formula  have  the  same 
asymptotics.  Both  approach  the  one  center  limit  7^^  or  7^^  as  Rab  0,  and 
approach  zero  as  Rab  — > 00.  We  decided  to  adopt  the  Ohno-Klopman  formula  to 
calculate  the  two-center  two-electron  integrals  and  the  two-center  core-electron  terms 

{^^A\-^\^A)  = (3.15) 

^AB 

where  7^^^^  is  obtained  from  Eqn.  (3.14),  because  it  has  been  shown  by  Zerner 
[37]  that  the  MN  formula  cannot  successfully  yield  both  singlet  and  triplet  excitation 
energies  for  aromatic  systems.  Adopting  the  OK  formula  may  need  higher  levels  of  Cl 
in  reproducing  singlet  spectra,  but  the  triplet  spectra  is  readily  reproduced.  Finally, 

{^^A\VB\^A)  = (3.16) 

where  is  the  electron  occupation  number  on  atomic  orbital 

Two-  and  three-center  nuclear-electron  integrals 

There  exist  two-center  nuclear-electron  integrals  {ij.a\Va\(Tb)  and  also  three- 
center  {^iA\Vc^A,B\cr  b)  ■ These  types  of  integrals  are  not  present  in  the  INDO  methods, 
because  the  only  off-diagonal  one-electron  integral  in  INDO  methods  are  the  reso- 
nance integrals  /?,  which  are  empirical  two-center  integrals.  By  using  the  product  of 
/3  and  the  overlap  alone,  the  contributions  from  these  two  types  of  integrals  are  not 
so  clear.  Only  in  an  average  manner  can  their  contributions  be  represented  by  ;5A. 
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In  the  XINDO  method,  these  two  types  of  integrals  are  evaluated  explicitly. 

For  any  integrals  which  are  multi-centered,  the  functional  form  assigned  to  them 
must  have  a proper  distance  dependence,  in  order  to  be  compatible  with  the  goal  of 
obtaining  good  geometries.  And  again,  from  the  fact  that  the  INDO/S  method  has 
to  use  different  weighting  factors  to  accompany  j3  in  order  to  obtain  good  spectra,  it 
may  be  better  to  explicitly  evaluate  the  different  parts  of  the  off-diagonal  one-electron 
matrix  elements. 

There  are  two  options  for  calculating  the  (i2a\Va  + Vb|(Tb).  The  first  option 
is  based  on  the  idea  that  no  extra  integrals  need  to  be  calculated,  instead,  they  are 
obtained  by  applying  the  Mulliken  approximation  [44]  using  the  diagonal  one-electron 
Hamiltonian  elements: 

(3.17) 

The  other  option  is  to  use  scaled  ab  initio  -H  Vb\ob)  integrals.  That  is, 

the  two-center  nuclear-electron  integrals  are  evaluated  from  first  principles  with  the 
orbital  exponents  from  Slater  rules  [29].  Then  from  the  one-center  nuclear-electron 
parameter  the  effective  exponents  are  evaluated.  Since  the  measure  of  electron- 
nuclear  attraction  is  the  orbital  exponent,  one  can  scale  the  first  principles  integrals 
by  the  ratios  of  the  effective  exponents  and  the  exponents  ^ from  Slater  rules.  In 
that  case  the  scaled  {^ia\Va  + Fb|ctb)  integrals  are  compatible  with  the  parameter 

%A- 

The  scaled  ab  initio  integrals  may  be  better  than  the  integrals  using  the  Mulliken 
approximation  when,  e.g.,  the  electronegativities  of  the  two  atoms  involved  are  very 
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different  since  the  Mulliken  approximation  will  force  the  charge  to  distribute  more 
equally  between  the  two  atoms.  Another  disadvantage  of  the  Mulliken  approximation 
is  that  it  does  not  have  the  right  asymptotics.  It  introduces  large  errors  when  the 
internuclear  distance  is  either  very  small  (where  it  may  not  matter)  or  very  large. 
Thus,  these  problems  may  be  avoided  by  choosing  the  scaled  ab  initio  integrals  and 
we  decided  to  use  such  integrals. 

The  three-center  integrals,  {f^A\ycJ^A,B\c^B),  are  important  for  obtaining  con- 
formational information.  These  integrals  can  become  quite  sizeable  with  increasing 
molecular  size,  thus  if  we  want  to  introduce  these  integrals  into  the  INDO  model, 
we  must  pay  attention  to  the  balance  of  the  model  by  introducing  three-center  two- 
electron  integrals  as  well.  With  the  intention  to  avoid  evaluating  extra  two-electron 
integrals,  we  will  adopt  the  Mulliken  approximation  for  these  three-center  integrals, 

(/m|  X]  ~ —^[{f-Al  X yclfJ'A)  + {(^b\  X ^c\<^b)\,  (3.18) 

where  (//.4I Vc|/^a)  and  (ctbII'cIo’s)  are  evaluated  according  to  Eqs.  (3.16)  and  (3.14). 
The  Mulliken  approximation  by  itself  is  not  rotational  invariant.  However,  Eqn. 
(3.18)  only  involves  empirical  atomic  parameters  7'  and  the  Ohno-Klopman  formula. 
Thus  there  is  no  angular  dependence  in  this  expression,  and  therefore  the  model 
remains  rotational  invariant. 

Using  the  Mulliken  approximation  can  introduce  errors,  as  mentioned  earlier, 
and  now  the  error  in  evaluating  the  three-center  integrals  when  the  third  center  (C)  is 
located  in  between  the  other  two  centers  (A  and  B)  is  especially  severe.  We  therefore 
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will  keep  this  in  mind,  and  we  shall  add  correction  parameters  to  lessen  such  errors. 
In  addition,  the  error  in  these  terms  should  cancel  against  the  similar  error  that  arises 
in  the  evaluation  of  the  two-electron  terms  as  discussed  below. 

Effective  core  potential 

Another  important  feature  of  the  ZDO  methods  is  that  they  are  all-valence- 
electron  methods  and  assume  that  the  core  electrons  are  very  localized  and  will  not 
participate  in  chemical  reactions  or  low  energy  electronic  spectroscopy.  Therefore 
the  core  electrons  are  not  included  in  the  calculations  other  than  modifying  the  core 
charges.  However,  the  core  electrons  also  repell  the  valence  electrons  which  prevents 
them  from  collapsing  into  the  core  region.  Hence  in  order  to  avoid  such  collapse  that 
also  would  conflict  with  the  Pauli  principle,  it  is  reasonable  to  include  a short  range 
repulsion  to  the  system.  This  may  improve  the  model  so  reasonable  geometries  may 
be  obtained  by  using  spectroscopic  parameters.  There  are  many  studies  related  to 
the  repulsion  between  core  and  valence  electrons,  and  the  use  of  the  effective  core 
potential  (ECP)  seems  capable  to  account  for  such  effects.  Eor  the  valence-electron- 
only  systems  the  following  form  of  the  ECP  was  derived  by  Zerner  [43]: 


{ECP) 


OtB^A 


XI  ^IXACtB^OtBOCB^t 
XI  ^fJ-AOlB^aB^OlBlJ-A- 


(3.19) 


In  the  last  step  in  Eqn.  (3.19),  F^gaB  is  approximated  with  orbital  energy  Sag  which 
may  be  obtained  from  ESCA  experiments  [54].  This  is  justified  from  the  fact  that 
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the  core  electrons  are  very  localized  so  Koopmans’  approximation  may  be  considered 
good  [43].  The  above  expression  shows  that  this  potential  is  proportional  to  the 
overlap  squared  and  hence  this  term  rises  rapidly  at  short  interatomic  distance  and 
becomes  negligibly  small  at  ranges  where  the  bonding  interaction  dominates.  How- 
ever, through  the  way  that  INDO/S  was  parameterized,  the  ECP  has  been  included 
in  the  diagonal  one-electron  integrals.  The  ECP  contribution  will  be  added  to  the 
off-diagonal  one-electron  matrix  elements  in  the  new  Hamiltonian  in  order  to  fulfill 
the  Pauli  exclusion  principle  and  the  pursuit  of  a versatile  semiempirical  method. 
In  addition,  the  inclusion  of  the  ECP  should  improve  the  tendency  of  semiempirical 
models  to  overbind. 

Off-diagonal  matrix  elements  of  the  kinetic  energy 

The  core  integral  (Eqn.  (2.62))  in  the  INDO/S  method  is  parametrized 

using  experimental  data  from  atomic  spectra.  Therefore  it  already  contains  the  ECP, 
nuclear-electron  attraction,  and  also  the  kinetic  energy.  We  do  not  need  to  put  the 
kinetic  energy  into  the  diagonal  part  of  the  one-electron  integrals  but  we  are  going  to 
include  it  in  the  off-diagonal  elements.  This  is  because  it  has  a very  different  distance 
dependence  than  does  the  overlap. 

Coulomb  integrals 

Although  we  construct  the  XINDO  method  over  STOs,  the  two-electron  in- 
tegrals are  calculated  by  using  the  LTO  density  matrix  and  the  LTO  two-electron 


integrals. 
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We  start  the  explanation  for  such  approximation  with  the  relationship  between 
LTO  and  STO: 

^ ^STO^-l/2^  ^3  20) 

and  therefore  the  LCAO-MO  coefficient  matrices  for  LTO  and  STO  will  be  different: 


= ^ x^™D,  (3.21) 

where  K is  the  size  of  basis  set,  and  C,  D are  the  molecular  orbital  coefficient 
matrices  over  STOs  and  LTOs,  respectively.  With  the  coefficient  matrix  and  the 
orbital  occupation  number  n one  can  construct  the  density  matrix: 

= CnC^  = A-i/2Ai/2CnC^Ai/2A-i/2 

= A“^/^DnD^A“^/^  = (3.22) 


i.e., 


-pLTO  _ ^l/2p5T0^1/2  23) 


The  density  matrix  is  then  used  for  calculating  the  Coulomb  integrals: 


aX 

o-A  ij 

= E Cf °E 

ij  (t\ 
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= (3.24) 

ij 

where  the  superscripts  s and  I stand  for  STO  and  LTO  respectively.  Indeed,  the 
Coulomb  integrals  in  Slater  type  orbitals  can  be  replaced  by  half-transformed  in- 
tegrals as  the  types  shown  in  Eqn.  3.24  together  with  the  Lowdin  density  matrix. 
Furthermore,  we  introduce  the  approximation 

~ (3.25) 

b ij 

~ [E  P^f^{iAf^A\3A^^A)  + E P^f^{^B^iA\iB^^A)]  (INDO).  (3.26) 

iAjA  is^A 

Based  on  the  fact  that  after  the  half-transformation,  the  sum  of  one-center  integrals 
tends  to  increase  and  the  sum  of  two-center  integrals  decreases,  leaving  the  total  sum 
approximately  unchanged  [55].  Eqn.  (3.25)  says  that  the  mixed  STO- LTO  integrals 
can  be  evaluated  over  LTO  alone,  as  is  traditionally  done  in  ZDO  methods.  The 
conclusion  that  the  two-electron  Hamiltonian  does  not  change  much  after  the  basis 
transformation  is  in  agreement  with  this  observation. 

In  the  INDO/S  parameterization,  these  integrals  are,  in  fact,  semiempirical 
Slater-Condon  parameters  [56]  and  745  as  in  Eqn.  (2.107).  In  XINDO,  the  one- 
center  Coulomb  matrix  elements  are  evaluated  the  same  way  as  they  are  done  in 
INDO/S,  except  the  Ohno-Klopman  formula  instead  of  the  Mataga-Nishimoto  for- 
mula is  used  for  7^5. 

The  NDDO  method  has  non-zero  two-center  Coulomb  terms  {31a(^bWa>^b),  and 
such  integrals  appear  to  be  beneficial  in  the  calculations  of  geometries.  Incorporating 
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these  integrals  can  be  done  without  the  efforts  of  evaluating  the  extra  two-center  two- 
electron  integrals.  By  applying  the  Mulliken  approximation  with  the  corresponding 
one-center  two-electron  integrals,  defined  as  the  following: 


= (3.27) 

A 

which  is  similar  to  Eqn.  (3.25)  but  only  spherically  symmetric  integrals  are  included 
in  and  this  types  of  integrals  will  be  used  to  construct  the  two-center  Coulomb 
integral  and  still  maintain  the  rotational  invariance. 

Then  the  two-center  Coulomb  integrals  can  be  calculated  according  to: 


J, 


t^ACTB 


i’^ilAliA  T '^O'biTb)- 


(3.28) 


Although  Eqn.  (3.28)  is  the  expression  for  the  two-center  two-electron  integrals,  it 
contains  not  only  the  two-center  terms.  If  one  observes  the  expression  of  which 

contains  contributions  from  two  types  of  integrals:  and  {^J-A(^B\|J'A(^B), 

where  center  B is  any  atomic  center  except  A.  Hence  ™ay  contain  three-  and 

four-center  contributions.  The  inclusion  of  such  integrals  is  consistent  with  the  treat- 
ment for  the  one-electron  part  of  the  off-diagonal  elements.  Errors  in  the  Mulliken 
approximation  here  should  cancel  against  similar  errors  in  the  one-electron  matrix. 

Exchange  integrals 

Due  to  the  ZDO  approximation,  the  only  exchange  integrals  for  which  the  half- 
transformation can  be  applied  are  where  both  (f)^  and  (j)i,  are  on  center  A. 
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Such  type  of  integrals  are  of  the  same  form  as  the  Coulomb  integrals.  However,  in 
general  the  half-transformation  cannot  be  applied  to  the  following  type  of  exchange 
integrals, 

= E (3.29) 

To  illustrate  this,  we  have  orbital  on  atom  A and  (j)^  on  atom  5,  and  if  we  apply 
the  half  transformation  from  STO  to  LTO  to  exchange  integrals, 

= E (3.30) 

uX  i/\ 

then  it  becomes  clear  what  the  problem  is;  only  one  of  the  two  orbitals  with  the  same 
electronic  coordinate  is  transformed,  and  one  can  not  apply  the  ZDO  approximation 
over  such  term  Therefore,  Eqn.  (3.30)  offers  no  further  simplihcation.  We, 

then,  keep  the  exchange  integrals  the  same  as  in  the  INDO  methods  without  further 
approximation  since  two-electron  integrals  change  little  between  STOs  and  LTOs 
[50].  It  is  possible  that  when  applying  the  Mulliken  approximation  to  the  two-center 
Coulomb  integrals,  at  the  same  time,  the  two-center  exchange  integrals  are  approxi- 
mated the  same  way  by  the  corresponding  one-center  terms,  and  with  an  assumption 
that  the  one-electron  correction  parameters  can  take  care  of  this  error  in  an  average 
manner. 

The  nuclear-nuclear  repulsion 

In  order  to  obtain  the  total  energy,  a nuclear-nuclear  repulsion  expression  com- 
patible with  the  electronic  part  of  the  Hamiltonian  is  needed.  The  choice  of  this 
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expression  also  has  the  influence  over  the  geometry  optimization  calculations.  In  our 
formalism,  the  asymptotics  of  core-electron  interaction  when  i?  — > 0 (one-center),  7 , 
is  much  smaller  than  the  Coulomb  repulsion  1/i?  — >■  00.  If  we  adopt  the  value  of 
l/R  for  the  nuclear-nuclear  repulaion,  the  repulsion  within  the  molecule  will  be  too 
large  in  comparison  with  the  attractive  interaction.  Therefore,  a natural  choice  is 
to  use  the  parameters  for  one-center  core-electron  attractions  accompanied  with  the 
Ohno-Klopman  formula  for  nuclear-nuclear  repulsion 


E, 


nucl—nucl 


E 

A>B 


7*  7* 


(( 


1 


(3,31) 


A comparison  between  the  theoretical  repulsion  integral  and  the  semiempirical 
value  is  shown  in  Fig.  (2.1). 

3.2.3  Details  of  the  Parameterization 

In  this  section,  the  integrals  that  are  used  in  XINDO  are  collected. 

Diagonal  one-electron 


+ J2^ab  + ECP,  (3.32) 

B^A 

where  is  the  atomic  core  integral  and  the  one  center  ECP.  The  atomic  core  inte- 
gral is  determined  from  the  average  experimental  energies  of  neutral  atomic  configu- 
rations. It  contains  the  one-center  kinetic  energy  and  the  one-center  electron-nuclear 
attraction  and  the  one-center  ECP  by  construction.  The  expressions  for  core  integrals 
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of  the  first  two  rows  of  atoms  are  given  by  Karlsson  and  Zerner  [57]: 

Uss  = -Is  for  H, 

for  Li, 

Uss  = -Is  - F°{ss)  for  Be, 

Uss  = —Is  — {Z*  — l)F°(ss)  — — - — ^-G^{sp)  for  B-F, 

o 

Upp  = -4  - (Z*  - l)F°(ss)  - for  B-F,  (3.33) 

where  F*"  and  G*  are  the  semiempirical  Slater-Condon  factors  [56],  and  Z*  is  the 
effective  core  charge,  and  U is  the  ionization  potential  of  the  electron  from  orbital 
0A:-  The  second  term  in  Eqn.  (3.32)  includes  the  attraction  between  nucleus  B and 
the  electron  on  A,  and  the  repulsion  between  the  core  electrons  on  B and  the  electron 
on  A.  In  XINDO  this  term  is  expressed  using  the  Ohno-Klopman  formula  Eqn.  (3.14) 
with  the  parameters  7^^  and  and  Eqn.  (3.16). 

Off-diagonal  one-electron 

H^aitb  — (M/i|  ~ — Va  — Vb  — ^2  Uc\ob) 

^ C^A,B 

— ~ {I^a\VaWb)  — (/^AlThlo'B) 

- E {HA\Vc\aB)  + (ECP)  + (3.34) 

C^A,B 
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where  is  a kinetic  energy  element,  evaluated  ab  initio.  The  integrals 

are  also  evaluated  ab  initio,  but  scaled  to  be  compatible  with  the  parameter  7^^ . The 

three-center  integrals  {tiA\Vc\oB)  are  evaluated  by  Eqn.  (3.18). 

The  last  term  in  Eqn.  (3.34),  is  a correction  term  which  is  introduced  to 

take  care  of  the  errors  arising  from  the  Mulliken  approximation,  unmodified  exchange 
integrals,  imperfect  STO  Hamiltonian,  etc..  The  Wolfsberg-Helmholtz  formula  is 
used: 

PK..,  = + 0K,),  (3.36) 

where  (3b^^  is  an  orbital  dependent  atomic  parameter. 

For  two  different  spin  orbitals  and  (j),,  both  located  on  center  A,  the  one- 
electron  integrals  contain  nothing  but  the  ECP  under  the  INDO  approximation: 


E 

(core) 


^IJ.AOCB^CtB  ^OB^A  ■ 


Diagonal  two-electron 

For  a closed-shell  system. 


(3.36) 


~ X]  — -(/i/lA^lCT^/i^)]  -I-  ^ PbV^IAB-  (3.37) 

B^A 

The  values  of  (fJ-A^A\tP^A(^A)  and  are  determined  by  the  corresponding 

semiempirical  Slater-Condon  factors  [56].  ^ab  is  evaluated  by  the  Ohno-Klopman 
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formula  Eqn.  (3.14)  and  INDO/S  one-center  Coulomb  integral  7^^: 

l^lA  ~ ^liA  ~ ^il-AI  (3.38) 


where  is  the  ionization  potential  and  is  the  electron  affinity  of  atomic  orbital 

^^lA  ■ 

Off-diagonal  two-electron 

In  this  catagory  there  are  two  types  of  two-electron  integrals,  and  Gfj,^^g. 


fiSTO 


f-iSTO 


■ ^MAb_(T  a-  T 

2 ^ '^o-bo'b; 


(3.39) 


(3.40) 


There  are  different  suggestions  [55]  concerning  the  form  of  the  overlap  elements  in  the 
Mulliken  approximation.  The  original  formula  by  Mulliken  [44],  which  is  the  orbital 
overlap,  is  used  in  evaluating  both  one-  and  two-electron  integrals  in  the  XINDO 
method. 

SCF  procedure  of  XINDO 


In  the  previous  section  we  have  introduced  all  the  terms  present  in  the  XINDO 
model  Hamiltonian.  For  a closed-shell  system  of  s-  and  p-orbitals. 


pSTO 
^ y- A fi  A 


U, 


y-A^A 


C'BjtA 


E 


9 p 

‘-’yAdCB^^OlB'^aByA 


°‘^(core)^^ 
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+ 

+ 

Bi^A 

(3.41) 

J?STO 

= 

“B(core)^'4 

— 

(3.42) 

pSTO 

= 

W. -{/'/.IC4K)- WlVsIai,)-  E %^«/'.4|V'cW 

+ 

{(7b\Vc\(7b))  Y S^^ac^acSaccTB^  ^ iiaHa  ^ ^obob 

“C(core)  ^ 

) 

— 

2'^^'.40-b7AB  + P^fiA<JB^ 

(3.43) 

Most  of  the  changes  are  made  for  the  off-diagonal  elements  of  the  Fock  matrix,  and 
this  is  clear  if  one  compares  Eqn.  (3.43)  with  Eqn.  (2.71).  Table  3.1  and  Table  3.2 
contain  the  Eock  matrix  elements  of  the  XINDO  and  INDO/S  methods. 
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Table  3.1.  Fock  matrix  elements  of  the  XINDO  methods,  in  (s,p)  basis  set 


F 

Two-electron 

^tJ'AfJ-A 

^<tb7;.^<xb+ECP 

PAAlif^A^AlfJ-A^A)  —^(^J-A^aI^AIJ^a)] 
+ Hb^a  PbbIab 

F 

^IJ'Al'A 

{^,^\-\V^-Va\va) 

0.0 

'LB^A{^J'A\VB\VA) 

ECP 

Two-electron 

\PnAyA[H^J-A^A\^A^A)  - {^iA^A\|J^A^A)] 

F 

^ ilA<TB 

{iia\  - IV^Ictb) 

ab  initio 

{^^A\  —Va  — Fb|cTb) 

ab  initio 

— Y.C^A,B{^^A\yc\(7B) 

-^Eac  «Ac[7;.^Ac  +7^bAc] 
+ECP 

Two-electron 

2 [’^tJ'AIJ-A  T J(Tb<^b]  2^^^a<^b1AB 

correction 
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Table  3.2.  Fock  matrix  elements  of  the  INDO/S  method,  in  (s,p)  basis  set 


Two-electron 

Y^Bi^A  Z*b1AB 

PAA[{^^A^A\^^A^A)  —\{^J‘A^A\^A^J‘A)\ 

+ J2b^a  PbbIab 

F 

^ il-Al'A 

{y,^\-\V^-VA\vA) 

0.0 

Y.Bi.A{l^A\yB\yA) 

0.0 

Two-electron 

lPnA‘'A[H^^A^^AWA^^A)  - {^J‘A^^A\^iAl^A)] 

F 

{^iA\  - 

{^J‘A\  — Va  — Vb\ob) 

^Wa+(3b] 

Two-electron 

~ 2^t‘AC'B'yAB 

Semiempirical  methods  of  the  ZDO  type  assume  an  identity  overlap  matrix  in 
the  process  of  diagonalizing  the  Fock  matrix  and  this  is  justified  by  considering  the 
Hamiltonian  transformed  to  the  orthogonalized  atomic  basis  set.  However,  traditional 
ZDO  methods  may  not  have  all  the  terms  for  the  orthogonalized  basis  set  included 
in  the  model  Hamiltonian  and  that  causes  difficulty  in  applying  the  ZDO  methods 
for  general  purposes.  Jug,  Thiel  and  many  others  have  tried  to  modify  the  ZDO 
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methods  by  adding  correction  terms  into  the  Hamiltonian  in  the  LTO  representation. 
The  XINDO  method,  on  the  other  hand,  constructs  the  Hamiltonian  in  the  STO 
representation,  similar  to  the  EHT  in  this  sense,  but,  of  course,  unlike  EHT,  our  new 
method  contains  not  only  one-electron  but  also  two-electron  integrals. 

In  XINDO,  the  Hamiltonian  is  built  over  a non-orthogonal  basis  set  and  yet 
keeps  the  evaluation  of  integrals  to  a minimum.  The  density  matrix  used  for  building 
the  two-electron  part  of  the  Fock  matrix  is  still  the  Lowdin  type  density  matrix  as 
illustrated  in  Eqn.  (3.24).  The  Roothaan  equation  (Eqn.  (2.18))  is  solved  with  non- 
identity overlap  matrix.  The  orbital  energies  are  obtained  from  the  solutions  of  the 
Roothaan  equation,  and  so  are  the  corresponding  molecular  orbitals,  and  thence  the 
density  matrix  which  corresponds  to  the  STO.  Total  energy  is  evaluated  based  on 
the  STO  density  as  the  Fock  matrix  is  created  over  the  STO  basis. 


£».a,  = 1 E + E 1 


(3.44) 


Here  the  nuclear-nuclear  repulsion  energy  is  evaluated  using  the  parameters  for 
electron-nuclear  attraction. 

In  such  an  arrangement,  we  include  the  overlap,  in  spirit,  into  the  XINDO 
Hamiltonian.  This  results  not  only  in  the  improvement  of  the  description  of  the 
antibonding  orbitals,  but  also,  through  the  diagonalization  procedure,  multi-center 
interactions  are  included.  Hence  this  treatment  may  benefit  the  calaulations  related 
to,  e.g.,  molecular  geometries,  multipole  interactions,  rotational  barriers,  etc.. 
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We  then,  back  transform  the  density  matrix  to  the  Lowdin  type  density  to  build 
the  two-electron  matrix  and  proceed  the  next  iteration  until  the  STO  density  from 
two  consecutive  SCF  cycles  differ  less  than  a given  tolerance. 

3.2.4  Similar  Work  in  the  Field 

Different  proposals  for  a multi-purposed  semiempirical  method  have  been  pre- 
sented by  others  groups  and  indeed  there  has  been  some  successful  attempts.  As 
mentioned  in  the  previous  section  most  of  the  them  tried  to  improve  the  ZDO  Hamil- 
tonian in  the  LTO  representation.  This  is  one  of  the  main  differences  between  XINDO 
and  the  other  methods. 

It  may  be  worthwhile  to  discuss  the  problems  that  possibly  exist  in  the  current 
ZDO  methods.  These  problems  are  related  to  the  effects  to  the  Hamiltonian  arising 
from  the  basis  set  orthogonalization  transformation.  As  mentioned,  it  has  been  shown 
numerically  by  Coffey  [51]  that  in  the  case  of  the  C2  molecule,  the  contributions  from 
the  effective  core  potential  (ECP)  together  with  the  orthogonalization  corrections, 
are  about  equal  in  magnitude  to  the  penetration  integrals,  ~ Vab-,  which  has 

been  stated  in  the  review  chapters  (see  Eqn.  (2.59),  (2.60),  and  (2.61)).  However, 
the  integrals  used  by  Coffey  were  all  theoretically  evaluated  using  STOs,  and  contain 
no  electron  correlation.  Such  integrals  are  ideal  for  geometry  optimization  but  not 
for  calculating  spectroscopic  properties.  There  is  no  reason  for  Coffey’s  argument 
to  hold  when  using  other  types  of  integrals,  e.g.,  the  two-electron  integrals  in  the 
INDO/S  method.  Furthermore,  a two-center  situation  is  far  easier  to  handle  than 


systems  with  three  or  more  atoms. 
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Many  other  studies  have  been  presented  concerning  the  orthogonalization  cor- 
rection [30,46,49-51,58-68].  The  recent  conclusion  from  these  studies  is  that  the 
orthogonalization  correction  should  be  included,  and  therefore  the  ECP  should  be 
taken  into  consideration  explicitly  as  well.  Jug  et  al.  has  presented  a semiempirical 
Hamiltonian  [30]  which  is  a model  at  the  INDO  level.  The  ECP  is  added  into  their 
one-electron  Hamiltonian,  and  the  orthogonalization  correction  term  is  obtained  by 
doing  a Taylor  expansion  of  the  overlap  matrix  (Eqn.  (3.12)),  and  collecting 

terms  up  to  the  second  order  in  S.  Similar  work  has  been  done  at  the  NDDO  level 
by  Thiel  et  al.  [50,67].  It  turns  out  that  both  the  ECP  and  the  orthogonalization 
correction  term  are  repulsive  so  that  the  consitituent  atomic  orbitals  are  raised  en- 
ergetically. However  they  are  still  subject  to  a symmetric  splitting  between  bonding 
and  antibonding  molecular  orbitals.  They  chose  the  parameters  properly  so  that  the 
energy  of  the  bonding  molecular  orbital  is  close  to  the  ionization  potential  (Koop- 
mans’  theorem),  which  also  result  in  an  energy  lift  for  the  antibonding  molecular 
orbital.  The  energy  of  the  antibonding  MO  has  been  underestimated  by  all  current 
semiempirical  methods.  It  has  also  been  reported  that  in  the  new  NDDO  model, 
three-center  interactions  are  included  from  terms  that  originate  from  the  overlap  ex- 
pansion. This  makes  the  prediction  of  molecular  conformations  somewhat  better. 
However,  there  are  still  no  calculations  of  electronic  spectra  reported  using  these 
models. 

The  differences  between  XINDO  and  the  latest  ZDO  models  by  Jug  [68]  at  the 
INDO  level  or  by  Thiel  [50]  at  the  NDDO  level  are  summerized  in  Table  3.3. 
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Table  3.3.  Comparison  between  the  XINDO  Hamiltonian  and  the  orthogonalization- 
corrected  Hamiltonian  by  Jug  et  al.  [68]  or  Thiel  et  al.  [50] 


Feature 

Orthogonalization 

corrected 

XINDO 

Basis  set 

LTO 

STO 

Multi-center  effect 

< 3 

> 4 

Penetration  integral 

from  expanding  the 
overlap  series 

from  assigning  7 
parameter 

Possible  flaw 

overlap  series  may 
not  converge 

parameter  needs  to 
be  further  improved 

Two-electron  integrals 

ab  initio 

data  from  fitting 
atomic  spectra 
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XINDO  ZDO  OAO  corrected 

Figure  3.1.  MO  schemes  for  XINDO,  ZDO,  and  orthogonalization  corrected  methods 


The  MO  schemes  in  Fig.  3.1  are  constructed  from  the  one-electron  matrix 
elements.  H^x  denotes  the  XINDO  elements  and  the  resulting  MOs  split  unsymmet- 
rically.  (which  is  similar  in  magnitude  with  H^x,  i-e.,  H^z  ~ 77n)  denotes  the 
one-electron  matrix  element  in  ZDO  methods,  and  stands  for  the  matrix  element 
which  is  corrected  with  the  atomic  orbital  orthogonalization  effect.  The  latter  two 
types  of  methods  yield  symmetrically  split  MOs,  but  the  repulsive  effect  from  the  or- 
bital orthogonalization  correction  raises  the  which  results  in  the  higher  energy 
for  the  antibonding  MO  compared  to  standard  ZDO.  The  reason  why  the  XINDO 
method  gives  an  unsymmetric  split  was  discussed  in  Eqn.  (3.2):  a result  of  including 
the  overlap  into  the  model. 

Prior  to  this  thesis  some  efforts  were  made  to  build  a new  Hamiltonian  within 
INDO/S  scheme  by  reformulating  the  functional  form  of  the  resonance  integrals  (5 
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[69,70]: 

= r"exp[-K,  (3.45) 

where  r is  the  interatomic  distance,  m,  a,  and  /3  are  parameters  in  the  model  to 
be  adjusted  to  obtain  a better  distance  dependence  of  the  resonance  integrals.  The 
equation  above  comes  from  the  understanding  that  the  distance  dependences  of 
and  are  different.  Since 


(3.46) 

and  according  to  Eqn.  (3.9), 

(3.47) 

Therefore  simple  forms  such  as  the  Wolfsberg-Helmholtz  formula  [23]  may  not  be 
sufficient  for  LTO  off-diagonal  one-electron  matrix  elements. 

This  model,  Eqn.  (3.45),  accurately  reproduces  electronic  excitations  because 
the  parameterization  scheme  is  exactly  the  same  as  in  the  INDO/S  method  with 
the  exception,  of  course,  of  the  resonance  integrals.  It  is  able  to  perform  molecular 
geometry  calculation  [70]  to  some  extent.  However,  multi-center  terms  are  still  not 
present,  thus  we  suspect  the  calculation  of  geometry  may  be  limited  in  accuracy. 


CHAPTER  4 

PARAMETERS  AND  OTHER  RESULTS 


4.1  Parameters 

All  the  parameters  in  this  method  are  chosen  as  orbital  dependent  atomic  pa- 
rameters. Two-center  bonding  parameters  are  avoided  because  it  is  a complicated 
work  to  select  a parameterization  set  with  the  capability  to  represent  all  sorts  of 
bonding  situations  for  different  molecules  that  contain  similar  atom  pairs. 

We  want  this  model  to  have  correct  asymptotic  limit  as  the  inter-nuclear  dis- 
tances become  large.  To  do  so,  we  may  again  consider  a homonuclear  diatomic 
molecule  with  bond  length  infinitely  long,  and  the  energy  expression  of  this  system 
is  simplified  to 
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(4.1) 


where  the  first  term  is  the  nuclear-nuclear  repulsion,  the  second  term  is  the  electron- 
nuclear  attraction,  and  the  last  term  is  the  electron-electron  repulsion.  {'y'y^)avg  is  the 
nuclear-electron  attraction  parameter  in  XINDO,  and  here  an  average  value  taken 
from  all  orbitals  in  the  atom.  Similarly  the  {'yA)avg  is  an  average  value  of  all  the  one- 
center  two-electron  integrals.  It  is  clear  from  this  energy  expression  that  at  infinite 
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distance,  the  binding  energy  becomes  zero.  We  also  want  the  potential  curve  of  this 
system  to  approach  zero  with  positive  slope  asympotically  {dEao/dR  > 0),  so  we 
differentiate  Eqn.  (4.1)  with  respect  to  the  inter-nuclear  distance  R: 
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from  which  we  obtain 
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so  we  may  deduce  the  following  relation 


{'lAjavg  ^ {'Ja)  avg^  (4.4) 

From  this  we  know  that  the  parameters  for  nuclear-electron  attraction  should  be 
larger  than  the  one-center  two-electron  integrals  in  average  to  give  the  XINDO 
method  a correct  asymptotic  behavior.  This  is  certainly  so  in  ab  initio  methods, 
for  example,  where  74  = 0.3633^  and  7^  = 0.5^  (^  is  the  exponent  of  a Slater  type 
orbital). 

Hence,  based  on  Eqn.  (4.4)  we  start  to  search  for  the  parameters  by  calculating 
the  homonuclear  diatomic  molecules:  H2,  C2,  N2,  O2.  Before  knowing  the  precise 
ranges  of  the  parameters,  we  do  XINDO  calculations  on  these  molecules  by  trial 
and  error.  For  each  of  the  diatomic  molecules,  we  first  guess  the  7^,  7^,  and 
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(3bp  parameters  for  each  atom  (except  H2,  where  only  s-parameters  are  applied). 
When  the  calculated  ionization  potentials  are  approximately  correct,  and  the  highest 
occupied  molecular  orbitals  have  the  correct  symmetries,  we  then  proceed  with  the 
geometry  calculations  with  the  same  parameters.  If  the  calculated  bond  length  is 
close  to  the  experimental  value  from  the  literature,  then  the  set  of  parameters  for 
this  atom  can  be  accepted  as  a good  initial  guess. 

The  next  step  is  to  tune  the  parameters  even  further  using  medium-sized  molecules. 
A set  of  molecules  are  used  for  geometry  optimizations  using  the  XINDO  method.  By 
the  aid  of  the  simplex  algorithm  [71-73]  we  can  minimize  the  deviation  between  the 
calculated  results  and  the  expected  values  and  finally  the  parameters  corresponding 
to  the  smallest  error  (locally)  are  obtained. 

Give  initial  guess  of  parameters 

^ XINDO  geometry  optimizations 

outputs  (geometry  and  orbital  energy)  read  by  perl  scripts 

Feed  the  outputs  to  the  goal 
function  by  perl  scripts 
to  evaluate  the  error  s 

i 

New  parameters  from  simplex 

Figure  4.1.  The  flowchart  of  the  procedure  of  obtaining  XINDO  parameters 
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Fig.  4.1  summerizes  the  prodecdure  that  is  used  to  obtain  the  C,  H,  N,  and  O 
parameters  for  the  XINDO  method.  The  more  molecules  included  in  the  training  set 
the  better  parameters  may  be  obtained.  Currently  the  training  set  contains  the  fol- 
lowing molecules:  methane,  ethane,  ethene,  ethyne,  benzene,  naphthalene,  pyridine, 
pyrazine,  pyrrole,  imidazole,  methylamine,  ammonia,  water,  phenol,  p-benzoquinone, 
furan,  propionic  acid,  ether,  methanol,  and  acetaldehyde. 

The  properties  that  we  use  as  the  criteria  for  XINDO  calculations  are  the  highest 
few  orbital  energies,  and  the  molecular  structures.  Good  parameters  can  have  the 
orbital  energies  close  to  the  experimental  ionization  potential  (Koopmans’  theorem), 
and  the  calculated  molecular  structures  similar  to  the  experimental  geometries. 

We  also  examine  the  parameters  for  different  sizes  of  molecules  to  ensure  that 
there  are  not  any  parts  in  the  Hamiltonian  that  introduce  errors  that  accumulate  as 
the  molecular  size  increases.  Such  “size  inconsistency”  would  cause  problems  when 
the  model  is  applied  to  cases  with  different  ranges  of  size  other  than  the  size  of 
molecules  in  the  training  set.  It  turns  out  that  no  problem  of  this  nature  exists,  and 
the  errors  among  various  sizes  of  hydrocarbons  are  found  similar. 

The  parameters  which  are  currently  used  in  the  XINDO  model  are  shown  in 
Table  4.1.  At  this  stage,  we  have  managed  to  obtain  the  parameters  for  H,  C,  N,  and 
O,  which  are  the  elements  that  constitute  most  organic  molecules.  The  parameters 
should  be  further  refined  and  therefore  they  are  called  “estimated”  for  the  time  being. 
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Table  4.1.  Estimated  parameters  (in  eV)  for  the  XINDO  model 


Atom 

(S4|v/i|s^)  =7^, 

{pa\Va\pa)  = i'aj. 

Pbs 

pbp 

H 

13.4170 

-37.6957 

C 

17.4366 

11.1128 

-42.3887 

-52.6009 

N 

18.2931 

11.61 

-51.6516 

-64.7970 

0 

17.6873 

12.7703 

-62.5875 

-97.8456 

The  parameters  of  one-center  two-electron  integrals  74  in  the  INDO/S  method 
are  shown  in  Table  4.2. 


Table  4.2.  7^  parameters  (in  eV)  in  the  INDO/S  method  (for  Ohno-Klopman  for- 
mula) 


Atom 

H 

C (both  s and  p) 

N (both  s and  p) 

0 (both  s and  p) 

lA 

12.85 

11.11 

12.01 

13.0 

It  is  clear  from  Table  4.1  and  Table  4.2  that  the  parameters  used  in  the  XINDO 


method  fulfill  the  relation  in  Eqn.  (4.4). 
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4.2  Test  Cases 

4.2.1  Geometry  Optimization 

To  check  the  performance  of  the  model  and  the  parameters,  several  examples 
are  calculated  and  illustrated  in  the  following.  First  of  all,  the  ability  to  optimize 
molecular  geometry  is  the  primary  goal  to  achieve  in  the  XINDO  model.  Therefore, 
some  geometry  optimization  outputs  are  reported  in  Table  4.3. 

The  calculated  results  seem  fairly  good  although  there  is  still  plenty  of  room  for 
improvement  as  said  earlier. 
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Table  4.3.  Optimized  geometries  with  the  XINDO  method 


Molecule 

Property 

This  work 

Experimental 

ethane 

C-C 

1.595  A 

1.536  A [74] 

C-H 

1.182  A 

1.102  A [74] 

ethene 

C=C 

1.316  A 

1.339  A [75] 

C-H 

1.145  A 

1.085  A [75] 

ethyne 

C=C 

1.132  A 

1.203  A [76] 

C-H 

0.997  A 

1.062  A [76] 

benzene 

C-C 

1.395  A 

1.397  A [77] 

C-H 

1.093  A 

1.084  A [77] 

naphthalene 

C-C 

1.392  A 

1.410  A [78] 

C-C 

1.408  A 

1.425  A [78] 

C-C 

1.360  A 

1.361  A [78] 

C-C 

1.429  A 

1.420  A [78] 

C-H 

1.065  A 

1.07  A [78] 

pyridine 

C-N 

1.328  A 

1.338  A [79] 

C-C 

1.398  A 

1.394  A [79] 

p-benzoquinone 

c=o 

1.246  A 

1.222  A [80] 

C-C 

1.469  A 

1.477  A [80] 

c=c 

1.347  A 

1.322  A [80] 

4.2.2 

We  examine  the  isomers  of  C3H3  to  demonstrate  the  conformal  properties  cal- 
culated by  XINDO.  The  INDO  methods  are  known  to  exaggerate  the  stabilization 
of  ring  conformers.  The  lack  of  three-center  integrals  might  be  reponsible  for  some 
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of  this  problem  since  in  that  case,  the  INDO  methods  cannot  properly  describe  such 
interactions. 

However,  the  ring  form  of  C3H3  is  anyway  predicted  to  be  more  stable  but 
energetically  it  is  only  1.08  eV  (25  Kcal/mol)  [81]  lower  than  the  linear  conformer. 
The  INDO  method  tends  to  overestimate  the  energy  difference  between  these  two 
conformers  because  of  the  tendency  of  the  INDO  methods  to  form  rings. 


Table  4.4.  Relative  energies  (in  eV)  of  CaH^  conformers  from  different  methods 


C3H3  conformer 

AMI 

INDO/S 

INDO/1 

XINDO 

Exp  [81] 

C2V  (Linear) 

0.0 

0.0 

0.0 

0.0 

0.0 

D3/J  (Ring) 

0.063 

-7.967 

-6.357 

-2.711 

-1.08 

Although  XINDO  still  overestimates  the  stabilization  of  the  D^h  conformer, 
there  is  an  obvious  improvement  over  the  INDO  methods. 

4.2.3  O. 

There  exist  two  conformers  of  the  ozone  molecule,  O3  at  the  Hartree-Fock  level 
of  theory,  bent  or  equallateral  triangle.  The  electronic  configuration  of  ozone  can 
be  described  by  three  valence  bond  structures,  consequently,  the  one-determinant 
description  of  the  SCF  level  calculations  is  not  able  to  give  the  correct  electronic 
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structure  with  only  one  minimum  at  117°.  Semiempirical  methods  bring  some  elec- 
tron correlation  into  systems  of  interest,  but  ozone  has  inherently  a two-configuration 
ground  state. 

The  multi-center  terms  included  in  the  XINDO  Hamiltonian  and  the  electron 
correlation  that  the  model  obtains  from  imitating  the  experimental  data  yield  an 
optimized  O3  structure  that  seems  very  good. 


Table  4.5.  Optimized  O3  geometry 


O3  conformer 

AMI 

INDO/1 

XINDO 

Exp  [82] 

0-0 

1.1604  A 

1.1438  A 

1.261  A 

1.278  A 

<000 

121.04° 

123.02° 

110.5° 

116.8° 

CHAPTER  5 

CONCLUDING  REMARK 


The  addition  of  penetration  integrals  (expression  (2.60))  has  made  it  possible  to 
obtain  good  molecular  geometries  by  using  the  semiempirical  two-electron  integrals, 
which  are  conventionally  used  only  for  calculating  spectra.  When  calculating  the 
geometries  without  the  ECP,  the  bond  lengths  tend  to  be  too  short  in  comparison  with 
the  experimental  values.  Such  results  are  expected  when  the  penetration  integrals 
are  applied,  as  pointed  out  by  Pople  [20].  However,  with  the  ECP  added  to  the  one- 
electron  part  of  the  Hamiltonian,  the  calculated  bond  lengths  are  much  improved. 

The  current  parameters  for  C and  O yield  the  dipole  moment  of  the  CO  molecule 
0.1712  Debye  pointing  towards  the  carbon  atom  (exp.  0.11  Debye  [83]).  This  is  an 
improvement  in  the  calculated  charge  distribution  compared  to  the  INDO/S  method. 
The  same  parameters  for  O also  yield  better  O3  geometry  compared  to  INDO/1  and 
AMI  methods  (see  previous  chapter).  The  geometry  optimization  of  the  staggered 
form  of  the  ethane  molecule  by  using  XINDO  gives  the  angle  109°  on  average  (110° 
and  108°).  The  multi-center  terms  that  are  included  in  both  off-diagonal  one-  and 
two-electron  integrals  results  in  a better  angular  description  and  also  the  charge  dis- 
tribution. Thus  from  the  promising  results  so  far,  this  method  seems  to  have  great 
potential.  Even  though  there  are  some  additional  terms  to  be  calculated  now,  most 
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of  them  are  one-electron  integrals,  and  the  integrals  are  evaluated  by  using  the 

corresponding  diagonal  elements  and  Jcjb<jb  according  to  Eqn.  (3.28).  There- 

fore, the  evaluation  of  these  new  terms  is  not  time  consuming,  compared  with  the 
time  needed  for  the  diagonalization  procedure,  which  is  the  same  as  in  the  INDO/S 
method. 

There  are  many  further  implementations  to  be  done  with  the  XINDO  model. 
First  of  all,  we  will  make  the  XINDO  method  able  to  perform  calculations  of  molecular 
spectroscopic  properties.  The  two-electron  integrals  that  are  used  in  the  XINDO 
method  are  the  same  as  in  the  INDO/S  method.  Therefore  we  expect  that  just  by 
finishing  the  compatible  Cl  program,  the  XINDO  method  should  be  able  to  give 
as  good  calculated  spectra  as  does  the  INDO/S  method.  Then  we  may  be  able  to 
calculate  molecular  geometries  of  excited  states  and  generate  excited  states  potential 
energy  surfaces. 

Up  to  now,  the  parameterization  of  the  XINDO  method  is  done  for  some  of  the 
first  row  elements.  It  is  far  from  a complete  model  if  it  cannot  be  used  to  calculate 
transition  metal  systems.  To  do  so,  we  need  to  formulate  the  Hamiltonian  with  proper 
expressions  and  parameters  for  a basis  set  including  s-,  p-,  d-  and  /-  functions. 

And  we  should  not  forget  that  the  XINDO  method  has  the  overlap  incorpo- 
rated. This  allows  us  to  develop  a semiempirical  method  which  is  parameterized 
over  an  extended  basis  set.  The  different  orbital  exponents  can  be  used  to  calculate 
different  sets  of  the  overlap  integrals.  By  applying  different  parameters  to  each  set  of 
exponents,  the  method  can  benefit  by  the  more  flexible  basis  sets.  Eventually  it  may 
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become  conceivable  to  calculate  Rydberg  state  electronic  excitations  using  extended 
basis  set  included  in  the  XINDO  method. 
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